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We consider tiie trial wavefunctions for the Fractional Quantum Hall Effect (FQHE) tliat are 
given by conformal blocks, and construct their associated edge excited states in full generality. 
The inner products between these edge states are computed in the thermodynamic limit, assuming 
generalized screening (i.e. short-range correlations only) inside the quantum Hall droplet, and 
using the language of boundary conformal field theory (boundary CFT). These inner products take 
universal values in this limit: they are equal to the corresponding inner products in the bulk 2d 
chiral CFT which underlies the trial wavefunction. This is a bulk/edge correspondence; it shows 
the equality between equal-time correlators along the edge and the correlators of the bulk CFT up 
to a Wick rotation. 

This approach is then used to analyze the entanglement spectrum (ES) of the ground state ob- 
tained with a bipartition ^UB in real-space. Starting from our universal result for inner products in 
the thermodynamic limit, we tackle corrections to scaling using standard field-theoretic and renor- 
malization group arguments. We prove that generalized screening implies that the entanglement 
Hamiltonian He = — log pA is isospectral to an operator that is local along the cut between A and 
B. We also show that a similar analysis can be carried out for particle partition. We discuss the 
close analogy between the formalism of trial wavefunctions given by conformal blocks and Tensor 
Product States, for which results analogous to ours have appeared recently. Finally, the edge theory 
and entanglement spectrum of p^ ± ipy paired superfluids are treated in a similar fashion in the 
appendix. 

PACS numbers: 74.20.Fg, 74.20. Rp, 03.67.Mn 



I. INTRODUCTION 

A. Motivation 

The Fractional Quantum Hall Effect (FQHE) is an 
archetype of strongly interacting many-body electronic 
systems. As the filling fraction is varied, various fully 
gapped phases of matter can be observed experimentally. 
These topological phases of matter exhibit spectacu- 
lar collective behavior, such as localized excitations with 
quantized fractional charge, abelian (and possibly non- 
abelian) fractional statistics, or protected gapless edge 
modes. Our understanding of the FQHE mostly relies 
on the "variational" approach pioneered by Laughlin [2] . 
His celebrated wavefunction is not the ground state of a 
physically realistic Hamiltonian (say, electrons in a mag- 
netic field with Coulomb repulsion). Yet, it describes 
accurately the FQHE at some specific filling fractions 
and it is commonly accepted that it belongs to the same 
topological phase as the real physical system. 

Since Laughlin's contribution, several types of trial 
wavefunctions have been proposed for the FQHE at var- 
ious filling fractions. These include for example hierar- 
chical states [3] , composite fermion wavefunctions j4j , or 
the important family of states given by conformal blocks 
introduced by Moore and Read (MR) [5] (states corre- 
sponding to Jack polynomials have also been proposed 
later [6] , they can actually be included in the latter fam- 
ily [?])• In this paper we focus on these MR trial wave- 



functions and their quantum entanglement properties. 

The idea that quantum entanglement [8i can help 
characterize different topological phases of matter has 
emerged over the past years. This approach has pro- 
vided valuable new insights, while traditional methods 
based on symmetry breaking and local order parameters 
are widely accepted to fail [T]. For instance, the topo- 
logical entanglement entropy of the ground state 
of a fully gapped Hamiltonian is one robust measure of 
quantum entanglement in a topological phase in two di- 
mensions. More precisely, a bipartition of the quantum 
system is defined when the Hilbert space factors into two 
parts % = Ha ® "Hb- For most physical systems, a natu- 
ral bipartition is given by a cut in position space in con- 
tinuous systems or between a row or a plane of sites in 
lattice models. The object of interest is then the reduced 
density matrix pA = tr-^^ (V'l- The von Neumann 
entropy of pA scales linearly with the area of the cut be- 
tween parts A and B |12j . and it contains a universal 
order 0(1) piece: the topological entanglement entropy. 
At the end of Ref. [TU], Kitaev and Preskill (KP) pointed 
out that the topological entanglement entropy arises nat- 
urally if one assumes that the reduced density matrix pA 
has the form of the thermal density matrix of a 1 4- 1- 
dimensional gapless chiral theory along the cut. 

The spectrum of the "entanglement Hamiltonian" 
He — — log PA is called the entanglement spectrum (ES). 
Li and Haldane (LH) studied the ES of quantum Hall 
systems numerically in Ref. |13j . and observed that it 
contains a low-lying part in which the multiplicities are 
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related in a universal way to the ones of the conformal 
field theory (CFT) describing the low-energy edge ex- 
citations. This low-lying universal part is usually well- 
separated from the rest of the ES by an entanglement 
gap. LH suggested that this low-lying part could be used 
as a diagnostic tool when comparing two ground state 
wavef unctions. In addition, they observed that for some 
specific trial states, such as the Laughlin wavefunction, 
the entanglement gap goes to infinity, leaving only the 
low-lying universal part. Although the work of LH re- 
lies on a bipartition that is a cut in momentum space 
(the so-called "orbital partition" [H]), which is not local 
[15j . their main observations remain valid for the more 
natural cut in real-space [TSHTT] . The real-space ES is 
expected to be the spectrum of a local field theory along 
the cut, in agreement with the insightful suggestion of 
KP: not only the multiplicities but also the eigenvalues 
of He = — log PA are the ones of a (perturbed) CFT 
along the cut, when the length of the cut is large com- 
pared to the mean particle spacing (which plays the role 
of an UV cutoff). The conjectured locality of the ES, also 
dubbed "scaling property" in [TS] , is discussed in greater 
detail in section II Dl 



The purpose of this paper is to provide an analytic 
framework that explains why these properties of the real- 
space ES hold for the large family of MR trial wave- 
functions. It involves a general construction of the space 
of edge excitations, and a precise analysis of the inner 
product in that Hilbert space. As a byproduct, we will 
arrive at the important conclusion that, assuming gen- 
eralized screening, i.e. short-range correlations only in 
the bulk (this will be discussed below), there is an iso- 
metric isomorphism (in the thermodynamic limit) be- 
tween the Hilbert space of the gapless edge excitations 
and the Hilbert space of the CFT used to construct the 
ground state trial wavefunction. This result is a precise 
"bulk/edge correspondence", which, stated loosely, says 
that the edge CFT and the bulk CFT are the same up to 
a Wick rotation. In particular, this rules out the possibil- 
ity of using non-unitary CFTs to construct FQHE trial 
wavef unctions, as those cannot be consistent with gen- 
eralized screening. Some version of this correspondence 
has long been expected [5] , although no precise statement 
about a general relation between the inner products in 
the space of edge states and those in the CFT has ever 
appeared in the literature. In Ref. [I8J, Wen provided 
an important argument that implies such a relation in 
the particular case of the Laughlin wavefunction, relying 
on the plasma mapping |2j ; the formalism we develop in 
this paper is strongly inspired by his. We will then ex- 
tend our ideas to tackle the real-space entanglement of 
the ground state. Let us emphasize that we will work 
only with wavef unctions, and do not address any ques- 
tion related to (physical) Hamiltonians in this paper. 

Previous steps towards an analytic understanding of 
the ES in quantum Hall systems include direct calcu- 
lations for the integer quantum Hall effect [Kl [H] or 
other free fermion systems pUH^ and rigorous results 



on the multiplicities for a large class of trial wavefunc- 
tions [23l [21] . We also note that an appealing physical 
scenario has been proposed in Ref. [5S]. It suggests to 
start from two pieces A and B of a topological phase 
which both support (counter-propagating) gapless chi- 
ral edge states, and then to glue them along the edge 
by switching on an interaction Hamiltonian that couples 
the systems A and B. Restricting the analysis to a cou- 
pling between the edges only, standard renormalization- 
group arguments yield the thermal form of the reduced 
density matrix suggested by KP. Despite its simplicity, 
this "cut-and-glue" argument relies neither on a specific 
wavefunction nor on a precise choice of the bipartition of 
the Hilbert space, and is therefore intrinsically different 
from the approach we adopt in this paper. 



B. Fractional quantum Hall effect and 
the lowest Landau level 

Let us first recall some standard notations for the 
many-body problem in the lowest Landau level (LLL). 
One considers N indistiguishable (spinless) charged par- 
ticles in a two dimensional surface pierced by a normal 
and uniform magnetic field. In this paper, this surface is 
either the plane or the sphere 5^ In both cases, we 
use complex coordinates to parametrize the surface: the 
plane is simply parametrized by z — a; + ij/ G C, while 
for the sphere of radius i?52, we use the stereographic 
coordinate z — 2Rs2 e^^'^ tan0/2, where {0,(l>) are the 
spherical polar coordinates. 

It is well-known that wavefunctions in the lowest Lan- 
dau level (LLL) correspond to analytic functions [5] . We 
therefore consider wavefunctions 



'^{zi, . . .,zn) 



(1.1) 



which are analytic in zi, . . . , zn, and satisfy the right 
statistics (either bosonic or fermionic) under the ex- 
change of the Xj's. On top of these requirements, 
\I'(zi, . . . , zn) must be normalizable: 



,zn)\ < (1-2) 



The measure e^^^'-^'-'d^Zi depends on the surface which 
we are considering. It can be computed directly by solv- 
ing the Landau problem for a single particle (see for 
the spherical case). The notation e^(^''^'^ comes from the 
fact that V{zi, Zi) can be viewed, in the plasma mapping 
[5] , as an electrostatic potential created by a ba ckgr ound 
charge (we will come back to this point in part HI I : 



-(A^* +2)log [l + |z|V(2i?s2)2] (sphere) 

V{z,z)= ■ 

(plane) 
(1-3) 

where £b is the magnetic length in the plane, and is 
the number of magnetic fluxes which pierce the sphere. 
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The Hilbert space corresponding to the LLL is finite 
dimensional in the sphere, but not in the plane. Indeed, 
for a single particle, a basis of normalizable analytic func- 
tions ^(z) is provided by the monomials z™, where m > 
can be any integer for the plane, while m S {0, 1, . . . , N,^} 
for the sphere. Keeping that remark in mind, in the rest 
of this paper, our notations allow us to treat the plane 
and the sphere on equal footing. 

The sphere and the plane both enjoy rotational invari- 
ance. The angular momentum in the plane is written Li . 
Note that z is the complex coordinate in the plane while 
z stands for a unit vector perpendicular to it. The basis 
of monomials z™ are eigenstates of 



(1.4) 



Meanwhile, with these conventions, the angular momen- 
tum if around the vertical axis of the sphere (also writ- 
ten z) acts on the monomials as 



L| . = - m) z"- 



(1.5) 



so the angular momentum on the sphere can always be 
related to the one in the plane. In particular, for N 
particles hf = - L^. 



C. Moore-Read construction 

We now recall how an Ansatz for the wavefunction 
^'(zi, . . . , zm) can be obtained by looking at conformal 
blocks in certain 2d chiral CFTs [5]. For a recent dis- 
cussion of the MR construction, and its implications for 
non-abelian statistics, see (26j. For basic CFT material, 
see E7H211. Let us consider the wavefunction 



N 



{N\Y[aiz,)\0) 



(1.6) 



where a(z) is a local operator (i.e. a primary field) in 
a given chiral CFT, acting on the vacuum |0). The 
charged vacuum (A''| will be defined precisely below; it 
carries a charge that is opposite to the total charge 
of n'^(-^i)j ill order for the correlator (1.61 to be non- 
zero. The chiral CFT and the field a(z) must be chosen 
such that ^'(zi, . . . , zjv) is single- valued: this implies that 
there must be a single fusion channel when one fuses the 
field a(z) with itself. The field a(z) can therefore be 
called abelian. Of course, for consistency, the function 
\I'(zi, . . . , z^v) must also be analytic and have the cor- 
rect statistics, which requires additional properties for 
the field a(z). For later convenience, the factor Zn is 
defined such that the wavefunction ^'(zi, . . . , Zj y) gi ven 



by the formula il.6l is normalized for the norm (1.2 1 



In the MR construction, the chiral CFT is chosen as 



the (tensor) product of two sectors CFT 



U{1) 



CFT 



the ^7(1) charge sector and an abelian field ipi^) in the 
statistics sector CFT^: 



a(z) 



X tp[z) 



(1.7) 



In this expression, ip{z) is a free chiral boson with prop- 
agator {(p{z)ip{w)) = — log(z — w), and normal order- 
ing is implicitly assumed in the exponential. The U{1) 
symmetry is generated by the transformations </?(z) i-> 
(p{z) + const. The correlator of the vertex operators 
{N\ n e*"^'^-)/ ^ |0) has to be invariant under these shifts, 
so the out vacuum (iV| must carry a f7(l)-charge propor- 
tional to N 



{N\ = lim 



1 



(1.8) 



This definition is standard in radial quantization of a 
CFT in the plane (see pTH^S] ). With this out vacuum 
(A^l, the correlator of the vertex operators is non-zero. It 

is equal to the Laughlin-Jastrow factor ni=i(-^« ~ ^jV^''^ 
leading to a trial state with filling fraction ly. The cor- 
relator in the statistics sector (nV'(^i)) depends on the 
choice of the CFT for the field ip{z). For example, the 
Laughlin wavefunction itself corresponds to the simplest 
case when ip{z) is the identity operator. Other possible 
choices of tpiz) include: a free fermion with propaga- 
tor l/{zi — Zj), leading to the MR (Pfaffian) wavefunc- 
tion ^ or minimal Fateev-Zamolodchikov parafermions 
[30] which give the Read-Rezayi (RR) series [31]. Other 
choices of the field ip(z) correspond to other wavefunc- 
tions which have appeared in the literature, for example 
the ones expressed in terms of Jack polynomials [5] . 

Like the U{1) charge sector, the statistics sector is as- 
sociated with some underlying symmetry, for example a 
Z2 symmetry in the case of a Majorana field generated 
by tp{z) —ip{z) (and more generally, a symmetry 
for parafermions). In that case, our definition of the out 
vacuum (1.8) must be completed by the insertion of a 



and the field a(z) is a product of a vertex operator in 



Z2-charge (Z^-charge) when N is odd [N ^ Omodfc), 
in order for the correlator (A^| Y[ o,{zi) |0) to be non-zero. 
The CFT for the statistics sector is always a rational one 
(there is a finite number of primary fields, which form a 
closed algebra under the operator product expansion). 



D. Entanglement of the trial states 

Because of their particular structure, the trial states 
given by conformal blocks have some very specific entan- 
glement properties, which are inherited from the under- 
lying CFT. In order to sketch some of these features, we 
first need to define the bipartition of the Hilbert space 
% = Ha (^H-b for which one computes the reduced den- 
sity matrix p^, or alternatively the Schmidt decomposi- 
tion 

I VI/)) =5]e-«-/2|*f))«|*f)), (1.9) 
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The ES depends on the bipartition and is, by definition, 
the set of pseudo-energies ^i's [13]. Throughout our pa- 
per we will use the notation with double rightangle for 
kets in the 2-1-1 physical Hilbert space, while simple righ- 
tangles are reserved for states in the CFT. The ground 
state l^*)) corresponds to the wavefunction (1.6). 



1. Real-Space Partition (RSP) 

A natural way to partition a system of itinerant parti- 
cles on some manifold M is to do a cut in position space. 
The single-particle Hilbert space "Hi is the L? space of all 
normalizable functions on M (in this section we do not 
require any analyticity condition). If we divide the man- 
ifold as M = MyiU A/s, then "Hi admits a decomposition 
as a direct sum "Hi = "Hia ffi ^is, where Hia {T~Lib) 
is the subspace of functions with support in Ma {Mb)- 
This simply means that any function /(r) can be written 
as /^(r) + /^(r), with /^(r) = (/^(r) = 0) if r ^ Ma 
(r (f. Mb)- This decomposition induces a corresponding 
bipartition of the A^-particlc space as 



N 
Na=0 



(1.10) 



where Na (Nb) is the number of particles in part A {B)- 
This bipartition is called real-space partition (RSP). 

In quantum Hall systems the RSP [121 US] is obtained 
by dividing the complex plane C where the coordinates 
Zi are defined (see section I B ) into two complemen- 
tary parts C = AU B. As is customary in the litera- 
ture, the bipartition of the plane C = AU B is usually 
taken such that the subsystem A is rotationally invari- 
ant (and simply-connected for simplicity): A is then a 
disc of radius R centered on the origin. The trial state 
I is usually an eigenstate of the angular momentum 
{Lf + L?) I ^E*)) = Li I 5')) , so the angular momentum of 
part A, noted Lf, is a good quantum number. The bi- 
partition of the Hilbert space with N particles thus 
takes the form 



N 
Na=0 



(1.11) 



with Na + Nb = N , and the different eigenvec- 
tors/eigenvalues in the Schmidt decomposition (5.2) can 



be classified according to the number of particles Na and 
to the angular momentum Lf . 

In general, for RSP there is a non-degenerate lowest 
pseudo-energy C(^AO:-^ix)) some values iV^o and L^q 
which depend on the system size N. We define A^, AA^^ 
and AL^ by subtracting off these values. 



2. Scaling property of the entanglement spectrum 

One can make a general conjecture, that is expected 
to hold for any ground state wavefunction in a 2 -I- Ic? 



topological phase, provided that the (connected) correla- 
tion functions of local operators evaluated in that ground 
state are all short-range: in each topological sector (i.e. 
for a fixed total anyonic charge in part A), the entan- 
glement Hamiltonian He = — log pA 'is isospectral to a 
"pseudo-Hamiltonian" that is local along the cut between 
A and B. Equivalently, in [15', this conjecture, dubbed 
scaling property, was stated as follows: for all AA^^ and 
AL^, as A^ — )• cx), the set of A^'s approach the energy 
levels of a "pseudo-Hamiltonian" that is the integral of a 
sum of local operators in a l-|-l(i theory defined along the 
cut between A and B. In particular, in phases of matter 
such as the FQHE states, there is a low-lying part (as 
observed by LH) that corresponds to a gapless sector in 
the 1 -t- Ic? theory. In general, the theory also has gapped 
excitations, which come with pseudo-energies larger than 
or equal to the entanglement gap. These gapped excita- 
tions are associated with Schmidt eigenstates that differ 
from the "cut ground state" {i.e. the Schmidt eigenstate 

^Wao )) ^^^^ smallest pseudo-energy), not 

only along the cut, but also far from the cut. Such states 
contribute with amplitudes that decay rapidly (exponen- 
tially) with the distance to the cut, hence the presence of 
an entanglement gap. The pseudo-energy levels that lie 
above the entanglement gap correspond to a mixture of 
gapped and gapless excitations. We note that for quan- 
tum Hall systems, this type of scenario has been dis- 
cussed in the case of the orbital partition in |32j . 

For some trial FQHE states, such as the Laughlin or 
MR (Pfaffian) states, the entanglement gap goes to infin- 
ity, leaving only the gapless low-energy part. One of the 
purpose of this paper is to explain why these wavefunc- 
tions, adn more genrally all the wavefunctions given by 
conformal blocks, exhibit this particular feature. Then 
we will also explain why, for these trial wavefunctions, 
the locality of the ES follows from the fact that all corre- 
lations are short-range, a property which is conveniently 
called generalized screening in this context, following 26|. 
We will find that, for the trial wavefunctions given by the 
MR construction, the ES is the spectrum of the Hamil- 
tonian of a perturbed 1 -I- Id CFT, and that this CFT 
is the one that underlies the ground state wavefunction. 
It is also the theory that describes the physical edge ex- 
citations, as we will show. In the simplest case where 
the CFT is perturbed only by irrelevant operators, the 
pseudo-energies are given by the eigenvalues of 



AC 



R 



Lo 



(1.12) 



where Lq = ALf -\-0{ANa) corresponds to the Virasoro 
generator of dilatations for the CFT in the plane. R is 
the radius of part A (for a circular region with perimeter 
2ttR) and v is some non-universal "velocity" . This veloc- 
ity has the dimension of a length, and is of the order of 
the mean particle spacing close to the cut, which is the 
natural UV cutoff in our problem. Thus the ratio v/R 
in our s caling is always of order I/^/Wa. The ellipses 
in (1.12 1 are terms of higher order in v/R ^ 1/\/Na, 
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which come from perturbing operators that are more ir- 
relevant. The precise dependance of the eigenvalues of 
Lq on AA^^ depends on the details of the CFT. For 
the Laughlin wavefunction at filling fraction v, one has 

Lo = ALf + More generally, the ES has to be 

discussed case-by-case, depending on what FQHE state 
one is dealing with. Depending on the CFT that un- 
derlies the state, some perturbing operators may or may 
not be present in the "pseudo-Hamiltonian" that gives 
the ES. It is also worth emphasizing that, just like for 
the physical energy spectrum in the presence of a real 
edge, the wavefunctions in the MR construction (with 
both the U{1) charge sector and a non-trivial statistics 
sector) usually have an ES with two branches of exci- 
tations (rather than one), with different velocities 
and Vjj; . More details about the perturbing operators are 
given in part |IV[ where we derive our main results on the 
ES. 

Another important aspect that we want to point out 
in this paper is the striking similarity between these par- 
ticular trial wavefunctions and other wavefunctions that 
are being used in condensed matter and in quantum in- 
formation: the Matrix and Tensor Product States (MPS 
and TPS). In that context, ground state entanglement 
properties have long been studied, with questions that 
are analogous to the ones that we tackle here. We will 
come back to this point in part |V] 



there is no reason to expect that the ES should be the 
spectrum of a local operator. Overall, if one looks at the 
whole spectrum, the PP is not local. However, for large 
Lf, we will see that the Schmidt eigenstates | "^f)) actu- 
ally correspond to the ground state configuration of Na 
particles localized in a circular cap centered on the north 
pole (and Ng particles in a cap centered on the south 
pole), and of edge excitations above this ground state. 
In that sense, and for Na/N 1/2 when N, Na — > oo, 
the part A {B) corresponds to the northern (southern) 
hemisphere. Thus, for these states, parts A and B may 
roughly be seen as spatial region, just like in the case of 
RSP. This observation can be made precise, and it has 
the consequence that the particle ES {i.e. the ES ob- 
tained from PP) at large Lf can be analyzed with the 
same techniques as for RSP. It leads to a similar result, 
namely that this part of the particle ES actually cor- 
responds to the spectrum of a local operator along the 
"cut" (the equator). Of course, as we highlighted in the 
last paragraph, the 5(9(3) rotational invariance of PP 
implies that the pseudo-energies depend on AL^ rather 
than on ALf, so this local operator needs to be different 



form ( 1.12 1, which is valid only for RSP. The particle ES 



will be discussed in greater detail in section IV D 



3. Particle Partition (PP) 

Another bipartition is natural for systems of identi- 
cal itinerant particles: the particle partition (PP) |33j . 
It is obtained by assigning a fictitious "pseudospin" de- 
gree of freedom for each particle. We label | ^)) or \B)) 
the two orthogonal pseudospin states. A spinless ground 
state is mapped into the larger space that includes pseu- 
dospin by assigning to each particle j the pseudospin 

state (^1 A)). + I B))^^ 7^2. The bipartition of this larger 

space is simply Ha <8) "Hb, where all the particles with 
pseudospin \A))- {\B))j) constitute part A (B). This 
definition is an extension of PP |15| , which includes dif- 
ferent particle numbers, unlike the usual one in the liter- 
ature [55] . 

In PP, rotational invariance is inherited from the in- 
variance of the full system. In particular, on the sphere 
the total angular momentum {L^)^ is also a good 
quantum number. In that case, the Schmidt eigenval- 
ues/eigenvectors can be organized in 50(3) multiplets. 
It can be shown easily that the Schmidt rank is the same 
in each Lf subsector both for RSP and PP [TSIIIS]. For 
PP, one can define Lf as the maximum total angular mo- 
mentum for the subsystem A, and A'^o the corresponding 
number of particles. Then we can define the quantum 
numbers AiV^ and AL^ (PP) by subtracting off these 
values. 

Unlike the RSP, the PP is not a local bipartition, so 



E. Structure of the paper 



Our paper is organized as follows. In section |II A[ we 
define some complementary notations for the trial wave- 



functions given by conformal blocks, and in section II B 
we argue that there is a natural and straightforward way 
of constructing the edge states which correspond to those. 
More precisely, we construct a linear mapping from the 
space of states in the CFT to the space of wavefunc- 
tions for the edge states. In part |III[ we give a detailed 
discussion of the screening property, and relate it to a 
boundary conformal field theory formalism. As a con- 
sequence, we derive one of the main result of this pa- 



per (section HID), which states that, assuming screen- 
ing, the quantum-mechanical inner products between the 
edge states are identical (in the thermodynamic limit) to 
the inner products in the conformal field theory. In other 
words, the linear mapping from section [II B| becomes an 
isometric isomorphism when the number of particles goes 
to infinity. We provide some numerical checks of this re- 
sult in section |IIIF[ Then, in part |IV[ we extend our ar- 
guments to the problem of computing the entanglement 
spectrum, and proving the scaling relation conjectured in 
\T5\ . Finally, in part |vj we emphasize the close relation 
between the wavefunctions given by conformal blocks and 
Tensor Product States, and discuss the relevance of our 
results in that context. 
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II. TRIAL WAVEFUNCTIONS FOR THE EDGE 
EXCITATIONS 

A. More structure for the CFT space 

1. The chiral algebra 

In general, there is a field a^{z) in the CFT such that 
the operator product expansion of with a is 



a{zi)a\z2) 



1 



(zi - Z2) 



2ha 



(2.1) 



Here we have introduced the conformal dimension ha of 
the field a{z). Since the field a{z) is chiral, ha is also 
its spin. In the quantum Hall literature, the conformal 
dimension ha is sometimes refered to as the spin per par- 
tide It is also related to the "shift" on the sphere 
[3] in a straightforward way: S = 2ha [23 ■ It is the sum 
of the conformal dimensions of the vertex operator and 
of the field (noted h^) 



h„. = 



2 



(2.2) 



The field a^{z) must carry a U{1) charge which is op- 
posite to the one of a(z), and similarly in the statistics 
sector. For example in the case of the Z^, parafermions, 
the field a{z) — e''^'^^^/^ (g) ^i{z) carries a charge 1 
(mod k), so must carry a charge fc — 1 (mod k). 
We have then 



(2.3) 



where 1p{zl)^lJ'^ {Z2) ~ l/(zi - ^2)^''* + . . . The fields a{z) 
and a^z) generate the chiral algebra A by the operator 
product expansion. 



2. Complex conjugation 



In this paper, we need to work with wavefunctions 



given by conformal correlators, such as (1.6 1, but also 
with their complex conjugate. For this purpose, it is use- 
ful to introduce an anti-chiral copy of the CFT, and an 
anti-chiral field 



a{z) = e'"^'^) X tPiz) 



(2.4) 



such that the complex conjugate of "^{zi, . . . , z^) is given 
by the following correlator in the anti-chiral CFT 

1 ^ _ 

[*(zi,...,zjv)]* = ^(iV|n«(^)|0) . (2.5) 

The field a(z), together with its conjugate a^{z), generate 
a copy of the same chiral algebra A. 












i ' a{z-i) 




\ • a{z'^ 


) \ 




i 

/ 

/ 

/ A 







(b) 



FIG. 1. The edge states wavefunctions are obtained (a) by in- 
serting contour integrals in the correlator as in formula (2.7 1, 
or equivalently (b) by computing the matrix element \2.d>\ 
with an out state {v\ in radial quantization. 



B. Edge excitations 



In this section we construct wavefunctions that resem- 



ble the "ground state" wavefunction (1.6), but which we 
interpret as the "edge excitations" . Note that we do not 
address physical Hamiltonians in this paper. Instead, 
the wavefunctions for the edge excitations are required 
to have the same short-range properties (cluster proper- 
ties) as the ground state (1.6), but they have different 
global properties, such as angular momenta. 

To construct these wavefunctions, we consider some set 
of fields (j)k{Ck) (fc = 1, . . . ,p) in the chiral algebra A, and 
we compute the correlator 



(2.6) 



i=l 



where the appropriate charged vacuum (iV'| is inserted 
(it must neutralize the total charge of the a(zi)'s and the 
(/)fc(C/c)'s). This correlator is a function of the Cfe's and 
of the Zi's. The short-distance properties as two or more 
of the Zi's come close to each other are inherited from 
the underlying operator product expansions. Therefore, 
they must be the same as the ones of the ground state 



wavefunction (1.6 1. This implies that the function (2.6) 
is analytic in the z^'s everywhere in the plane, except 
possibly at the points Cfc's, where it can have some sin- 
gularities. The function is always single valued, and the 
possible singularities at the C^'s are poles. In other words, 
the correlator (2.6) is meromorphic, and it is not yet a 



valid wavefunction for particles in the LLL. However, let 
us consider instead the contour integrals 



{N 



•^Cp -^Ci ^_-|^ 

(2.7) 

where the contours Ci , . . . , Cp encircle all the z^ s as 
shown in Fig. [Tja (and the contours are radially ordered: 
Ci is encircled by C2, and so on), and mk G Z. In cer- 
tain cases, the correlator ( |2.7| ) can be zero (for instance, 
this may happen if some of the m^'s are negative, when 
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the correlator (2.6) is analytic at Cfc oo). The con- 
tours can now be deformed, without changing the value 
of expression ( |2.7| ). In particular, they can be taken as 
circles with arbitrary large radii (i.e. they can be "de- 
formed around infinity"). Thus the expression ( |2.7[ ) is 
analytic in the Zj's in the whole complex plane, just like 
the ground state ( 1.6 ). If it is non-zero, it can be used as 



a FQHE trial wavefunction. Another advantage of using 
the contour integrals (2.7) rather than the meromorphic 



functions (2.6) is that one gets angular momentum eigen- 



states, which will be more convenient for the manipula- 
tions below. 

If we introduce the Hilbert space of the CFT 'Hcft, 
which is a (irreducible) module over the chiral algebra A, 
then (2.7) can be reformulated as follows. To each state 
we associate its dual (v\ S ^cft- Then we 



\v) e HcFT. 

construct the correlator 



AT 



(2., 



which we use as a wavefunction when it is non-zero (these 
wavefunctions appeared previously in [3U [33] ) . Thus, by 
definition, we have a linear mapping from the (dual) CFT 
Hilbert space T-Lq-pt the space of edge states. This 
mapping is in general not injective. 

We now consider two concrete examples, the Laugh- 
lin wavefunction and the MR (Pfaffian) wavefunction, to 
illustrate how this construction works in practice. 



1. First example: edge states for 
the Laughlin wavefunction 

For the Laughlin wavefunction, we have a{z) = 
gii^(z)/Vi^_ There is no statistics sector. The chi- 
ral algebra A is generated by the vertex operators 
giv(z)/v^ and e-''^(^)/'^ and by the operator prod- 
uct expansions. In particular, the U{1) current 
idifiz) is generated by e^'^^^i^/^e^^'^t^^)/^ - (zi - 
Z2)~'^^'' [I + (zi - Z2)idip{z2)/^/l^ + . . .]. The modes of 
the U{1) current J„ = ; 
mutation relations 



' d( ^^^idip{() satisfy the com- 



(2.9) 



As claimed in section [TQ the ground state wavefunction 
(1.6) is, up to the normalization factor. 



N 



(JV| J|e^^fe)/v/^|0) = Y[{z,~z, 

i=l i<j 



(2.10) 



The neutral edge states are obtained by exciting the out 
vacuum (A^| with the positive modes J„ {n> 0): 



N 



N 



N 



N N 



N 

TV N N 

^ E E E n(-^ - -.o'^^ (2-11) 

fc=i p=i 1=1 i<j 

and so on. The negative modes J„ [n < 0) do not 
need to be considered, since they annihilate the out vac- 
uum (7V|. In general, the positive mode J„ of the U{1) 
current J(z) produces the corresponding sum of pow- 
ers -^J2zi' EH ■ These wavefunctions are the well- 
known (neutral) edge states for the Laughlin wavefunc- 
tion |18j . Equivalently, one could have obtained the same 
space of edge states (up to some change of basis) by 
acting with the modes = dC C"-i+E7e*v'(C)/v^ 
and (at)„ = § dC ("-i+ 277 e-*^(?)/v^. The advan- 
tage of working with the modes of the C/(l) current 
here is that one has a canonical basis of the CFT 
states with the U{1) charge corresponding to N parti- 
cles: {\N) , J_i \N) , J_2 \N) , ,...}. It would be 
more painful to write such a basis in terms of states of 
the form a_„j . . . a_np{a'^)-mi ■ ■ ■ (a^)-mp \N}, although 
in principle nothing prevents us from doing that. 

Finally, although the discussion in this section is about 
neutral excitations (the number of particles N is identical 
to the one in the ground state), it is not difficult to extend 
it to the case of charged excitations. To get those, one 
needs to add some particles. For example, for a single 
particle added, the correlator (-/V + 1| nil~i^ ^(-^0 |0) is 
simply the Laughlin wavefunction with iV + 1 particles. 



2. Second example: edge states for the Moore-Read 
(Pfaffian) state 

The second example we consider is the MR (Pfaffian) 
wavefunction, which corresponds to a{z) — e^'c'^'^^l ^ x'l}}^ 
where is a free (Majorana) fcrmion field with propaga- 
tor {tp{z)^{'w)) = l/(z — w). For even particle number 
N , the ground state is then, up to normalization. 



AT 



N 



(ivin«(^.)io) - Wiz^-z^f" pf 



i=l 



Zk - Zl 



l<k,l<N^ 

(2.12) 
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In addition to the J„'s appearing in the U{1) sector as in 
the LaughUn case, the chiral algebra contains fermionic 
modes tpn = ^ f d( ^"^2 'ip{(), where 71 G Z + i, with 



(2.13) 



The neutral excitations obtained from the U{1) sector 
generate edge states which are similar to the ones of the 
Laughlin wavefunction. In the Majorana sector we have 
instead (for N even) 



was used later to show that the quasi-particles also obey 
(abelian) fractional statistics under adiabatic exchange. 
Wen used the screening property, coupled to an electro- 
static argument (the method of images |37] ) to construct 
the edge theory of the Laughlin states from a microscopic 
point of view |18j . 

In this part, we use a generalization of the screening 
property of Laughlin's plasma, discussed recently in full 
details in This "generalized screening assumption" 
is formulated as follows. The normalization factor of the 



(iV| VriiV'^ 



Pf 



N 
J = l 

( 





|0) = Y[('.-Zj) 



(2. 







wavefunction (1.6 1 can be viewed as the partition func- 
14);ion of a two-dimensional system of itinerant particles 

subject to some interactions and in a background elec- 
\ trostatic potential V{z,z) 



-1 



21-22 




'-N 







l<j 



"1-5 "2-5 
sgna X ^,(1)^^,(2)' 



""'(f - 1)! aTs^ (^'^(3) - ^.(4)) • • • (^a(Ar-l) 



where the Pfaffian of the (TV + 2) x (iV + 2) 
skew-symmatric matrix is the free-fermion correlator 
(V'ni'0n2'0('2^i) ■ • ■ V'(-2Ar)) that must be evaluated using 
Wick's theorem. Of course, more fermion modes can be 
inserted in the correlator to generate other edge states, 
and one recovers the wavefunctions constructed in |36j . 
The MR wavefunction with odd particle number also fits 
naturally in that framework, by inserting one fermion 
mode V'1/2 in the out vacuum (A^|, as discussed in sec- 
tion ITCl 



III. SCREENING 

It is well known that the normalization factor Z^q of the 
Laughlin wavefunction is exactly equal to the partition 
function of a two-dimensional one-component plasma in 
a background potential V{z,z) 



Zn — 



m 



N 



■^Z- f^i V{zj.Zj)+u-^ Y.k<l log |zfc-2,|' 



(3.1) 

This exact relation, usually refered to as "plasma map- 
ping", has been exploited in various ways in the litera- 
ture. A key point in the plasma mapping is the obser- 
vation that for an inverse filling fraction lower than 
about 70, the plasma is in a screening phase. This prop- 
erty was highlighted already by Laughlin, who used it 
to derive the fractional charge of the quasi-particles. It 



Zn — 



1 

M 



N 

n- 



y(z^ 



U^z., 



N 



{N\Y[aiz,)\0) 



(3.2) 



Contrary to the case of the Laughlin wavefunction, in 
general this partition function is not the one of a Coulom- 
bic plasma (i.e. involving only Coulomb two-body inter- 
actions). Instead, it is argued in [35] that the partition 
function (3.2) should in general be viewed as the one of 



a perturbed CFT (in a grand-canonical description, as 
we do in section III B I . Then, two situations may oc- 



cur, depending on the IR fixed point towards which the 
perturbed CFT is sent under the RG flow: either (i) 
the IR theory is massive, that is, all the connected cor- 
relations of local fields decay exponentially, or (ii) the 
IR theory contains massless modes and therefore long- 
range (power-law decaying) connected correlations. We 
say that "generalized screening" holds if the situation (i) 
occurs. It generalizes the case of the screening phase for 
the Laughlin wavefunction, which contains only exponen- 
tially decaying connected correlations. 

In general, there is no known analytical argument that 
allows to discriminate between the situations (i) and (ii) . 
Instead, one usually has to rely on indirect numerical 
checks of some of the consequences of the generalized 
screening assumption (i) . It is plausible that the property 
(i) holds for the MR (Pfaffian) wavefunction, and maybe 
for other states like the k — 3 RR state. Also, some 
general arguments have been given in |26j which show 
that the generalized screening cannot hold in some cases, 
as it would lead to contradictions (in particular in the 
case of non- unitary CFTs). 

In what follows, the "generalized screening" property 
(i) is assumed to hold, and we explore some of its conse- 
quences. This part of our paper is devoted to reformu- 
lating the screening property in the language of bound- 
ary CFT, and to using it to make a precise statement 
of the long expected "bulk/edge correspondence". The 
arguments in this section may be viewed as the natural 
generalization of Wen's microscopic theory of the edge 
excitations, which in its original formulation was appli- 
cable only to the Laughlin state [IS]. 
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FIG. 2. In the thermodynamic limit (A'^ — >■ oo), the particles 
fill a circular droplet of radius R. For the background poten- 
tial ( |1.3[ l corresponding to the plane, the density is uniform 
(a), for the sphere it is not uniform after the stereographic 
projection (b). 



Finally, for technical reasons, in the next sections we 
will often have to refer to the radius R of the circu- 



lar droplet for the rotationally invariant potentials (1.3 1 



The radius is fixed by the requirement that the number 
of particles inside the droplet is 



TV = 



Po{z,z)(fz. 



(3.4) 



\z\<R 



In the case of the constant background charge for the 
plane (1.3), this of course leads to the well-known radius 
R — y/2N/i'£B. On the sphere, the particles fill some 
spherical cap that is mapped onto the droplet by the 
stereographic projection, and the radius R depends both 
on the radius of the sphere and on the size of the cap. 



A. The droplet 

Consider the distribution of particles corresponding 
to the partition function (3.2 1 given by the normaliza- 
tion factor of the ground state wavefunction. In the 
large limit, the particles fill some domain in the plane 
C (where the coordinates z^'s are defined) called the 
"droplet". When the potential V{z,z) is rotationally 
invariant, as in (1.3 1, the droplet is circular and cen- 
tered on the origin (see Fig. [2]). At large distances (that 
is, much larger than the mean particle spacing, which 
is of order of £b in the plane) the average density of 
particles is well approximated by a continuous function 
(5(^^(2; — Zi) — >• po{z,z). It can be shown easily, for 
instance with a saddle-point approximation (a fully de- 
tailed derivation can be found in |3H])j that the density of 
particles is fixed by the U{1) charge sector only, and that 
it is equal (at scales >• £b) to the "background charge 
density" 



Pq{z,z) 



AF(z,z) 



(3.3) 



where A = Adzdg is the Laplacian. This background 
charge density is constant for the quadratic potential 
V{z,z) corresponding to the plane in formula (1.3). It 



is not constant for the sphere because the stereographic 
projection does not preserve the volu me (F ig. [2]). Let 
us emphasize the fact that the relation (3.3) has nothing 



to do with the screening property described in the pre- 
vious section. In particular, it holds in the crystallized 
phase for the Laughlin wavefunction, as well as in the 
screening phase, as long as the variations of the "back- 
ground charge" —^AV{z,z) occur on distances much 
larger than the mean particle spacing, such that a contin- 
uous description of the particle density is meaningful. In 
this paper, we will always be in the latter regime, where 
the the background potential V{z, z) varies slowly. For 
the case of the plane in (1.3), this is obviously true since 



the background charge does not vary at all, and for the 
sphere it varies on scales of order ^ Rg^ (the radius of 
the sphere) while the mean particle spacing is of order 
~ i?52 /VN, so their ratio vanishes when N ^ oo. 



B. Screening as a conformal boundary condition 

In this section we assume that the generalized screen- 
ing property holds, and we interpret its consequences at 
the edge of the droplet using the language of boundary 
conformal field theory (for a classical discussion of bound- 
ary CFT, see [39l HQ]). We will have to make a certain 
number of technical choices in order to state our argu- 
ments. The technicalities, however, should not prevent 
the reader from catching the basic idea, which is simple. 
Let us summarize it here. We start from a non-chiral 
massless theory CFT (g) CFT defined on a 2d surface. A 
subset of this surface, the "droplet" , is filled with parti- 
cles. In the field theory, these particles are equivalent to 
a perturbation of the form J a{z)a{z)(P z that is turned 
on inside the droplet (but not outside). This new term 
in the action drives the perturbed theory to a massive 
IR fixed point. That is the screening assumption. Now, 
outside the droplet, there is no perturbation, so the field 
theory CFT (X) CFT is still massless. Inside te droplet, 
all the correlation functions decay exponentially, and the 
correlation length is zero at the IR fixed point. Thus, we 
are left with a non-chiral CFT outside the droplet, and 
the local fields in this theory must satisfy a local boundary 
condition along the boundary of the droplet. The pur- 
pose of this section is to find this boundary condition. 
Its consequences for the edge theory of the quantum Hall 
states will be discussed later. Now let us turn to a more 
detailed formulation of this argument. 

It is more convenient to work in the grand-canonical 
ensemble as in [26] . Also, to avoid phase factors and nor- 
malization constants that would obscure the argument, 
we work on the cylinder C (see Fig. [3| parametrized by 
w ~ X + iy where y is identified with y + L. We imag- 
ine that the left half-cylinder Ci (Rew < 0) is filled by 
the particles in a uniform neutralizing background with, 
say, constant density pq. The way to treat this back- 
ground charge was discussed in [5]. Eventually, one can 
regularize this at a; — ^ — oo by integrating over the region 
—A < Re w < 0, such that the total background charge is 
finite, and then take A — oo in the end of the calculation. 
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The partition function of this system is 

(3.5) 

where the first exponential generates integrals over C; 
of pairs a{w)a(w). Following |26j . we look at this term 
as a perturbation of the action of the non-chiral theory 
CFT(g)CFT by the local operator a{w)a{w) in the region 
Ci. The coefficient A is included such that the term in 
the exponential is dimensionless. It can be tuned arbi- 
trarily, giving different weights to the terms with differ- 
ent number of pairs a{w)a{w). This will be important 
later. For now, note that because of the charge neu- 
trality of the correlator (3.5), only one term (the one 



with the U{1) charge that is exactly opposite to the to- 
tal background charge contained in region Ci) in the ex- 
pansion of the exponential actually contributes to the 
partition function Z(pq). The correlator is evaluated on 
the cylinder rather than in the plane, hence the sub- 
script C, and the propagator of the free chiral boson is 



Mwi)ip{w2))c = - log 



tanh 



The argument 



will not depend on the details of the background charge 
though, the important point here is simply that there is 
a non-zero density of particles po in the left half-cylinder 
Ci. 

Imagine that we want to compute the correlation func- 
tion of two operators a^wi), a'''{w2) in the right half- 
cylinder Rex > 0, as shown in Fig. [Sj 



a^wi)a^ {W2) e 



^fc, d^w a(w)a(w) ^-i^ f^^ d^w'[ip{w')+<p{w')] 



ZiPo) 

(3.6) 

(The reason why such correlators of fields outside the 
droplet are interesting will become clear below.) The two 
operators a^{wi) = p'Mwi)/V^ ^ i/jHwi) and a^w^) = 
'^('^)/V^ ^^\w:^) add a C/(l) charge to the correlator 
The latter is still non-zero though, because there 



(3.6) 



is again a term in the generating function of the integrals 
of pairs a{w)a{w) that has exactly the right U{1) charge 
required to ensure the global neutrality of the correla- 
tor. Similarly, in the statistics sector, the global neutral- 
ity (under transformations for example, if the statis- 
tics sector is generated by a Z/j-parafermionic current) is 
forced by the exponential that generates the pairs. 

To get some insight, let us first consider the C/(l) 
charge sector only. We use an electrostatic language, 
assuming screening in the "plasma" that fills the left 
half-cylinder C;. If the left-right symmetric operators 

"electric charges", then the operators e'^^v(.'u') I ^/^ and 
g±zip(u))/ themselves contain both electric and magnetic 
charge. In the plasma, the magnetic charge is confined, 
that is correlators of operators carrying magnetic charge 
fall exponentially (and their expectations in an infinite 
system vanish, so there is no need to subtract their dis- 
connected parts). The electric charge is screened in the 




• * • * 



\ y ' » » " • _ •* 



(a) 




(b) 




a<{w2) 



(c) 







FIG. 3. (color online) (a) Particles in a neutralizing back- 
groung on the left half-cylinder. We are interested in the 
correlation function of (wi) and a^{w2) in the right part, 
in the presence of all the particles, (b) Assuming screening, 
and in the hmit when the density of mobile charges in the left 
half-cylinder goes to infinity, an operator a^(wi) brought to 
the boundary from the outside is equivalent to an operator 
■ aljUi). (c) In this limit, the correlation function of a'{wi) 
and ci' {W2) in the presence of all the mobile charges in the 
left is equivalent to the two-point correlation function in the 
right part, with the boundary condition a{x,y) = a^(a;,y) at 
a; = 0. 



plasma, so correlators of electric charges fall also expo- 
nentially, with a correlation length of order ~ 1/ ^/po set 
by the density. This then implies that, when we take the 
operator a^{wi) to the boundary of the plasma from out- 
side (Rewi — >■ 0), it has its electric charge neutrahzed, 
leaving its magnetic charge. Thus, at the boundary of the 
plasma (Re wi =0), when inserted in the correlator (3.6 1, 
the operator e-*'^("'i)/v^ can be replaced by e'^^^^/^'" 
when the density po goes to infinity. Although this re- 
placement apparently violates charge neutrality because 
e-M^i)/V^ and e"^^^^/^ have opposite electric charge 
(but the same magnetic charge), it is valid inside the cor- 
relator (3.6), once again thanks to the exponential gen- 



erating integrals of pairs a{w)a{w) that ensures global 
charge neutrality. Thus, for the U{1) sector, we have the 
boundary condition along the imaginary axis (Re w = 0) 



-i^{w)j sjv 



(3.7) 



Now let us come back to the case of the full opera- 
tor a^(wi) = e~*'^(™i^/^ (g) V'^(wi)- Again, one brings 
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(wi) to the boundary (Rewi — )■ 0) from the outside. cyUnder 



In the correlator (3.6), it can fuse with one of the fields 
a{w)a{w), leaving only the field a{w). Assuming expo- 
nentially decaying correlations inside Ci, the "pairing" 
between a^{wi) and a{'w) can occur only on distances 
^ I/^/Pq. Therefore, when the density po goes to infin- 



ity, we are left with a CFT in the region Re it; > 0, where 
the fields are constrained by the boundary condition 



a\w) — a{w) 



(3.8) 



along the imaginary axis Re w — 0. This is a general- 



ization of the boundary condition (3.71 that includes the 
statistics sector. For example, when the statistics sec- 
tor is generated by a Z^-parafermion current, the con- 
straint (3.8) is a boundary condition on the Z^-current 



that was discussed in [HI . Again, the boundary con- 
dition (|3.8[) apparently violates charge neutrality, but the 



correlator (3.6) is still globally neutral thanks to the gen- 



erating function of integrals of pairs a{w)a{w). Strictly 
speaking, the replacement a^(?i;i) by a(wi) close to the 
boundary is only correct up to a multiplicative constant, 
which depends on the value of A in (3.6). Such a mul- 



tiplicative constant should also appear in the boundary 
condition (3.8). However, the coefficient A can always be 



tuned such that the multiplicative constant is 1, leading 
the simplest form of the boundary condition (3.8). 



The calculation of the correlator (3.6) then boils down 



to the one of the two-point correlator 



(3.9) 



in the domain Cr, with the boundary condition (3 



This is a considerable simplification of the problem. We 
will use this trick again in section [lII E| to compute equal- 
time correlators at the edge of a quantum Hall system. 



The boundary condition (3.8) for the non-chiral CFT 



outside the droplet is the main result of this section. It 
will play a crucial role in the rest of this paper. We 
obtained it from the specific form of the perturbation 
of the action / a{w)a{w)cPw inside the droplet, and as- 
suming that the generalized screening hypothesis holds 
(exponentially decaying connected correlations inside the 
droplet). The boundary condition (3.8) is a local con- 



straint along the boundary. It is a conformally invariant 
boundary condition: it is invariant under conformal map- 
pings w I— > f{w) of the domain outside the droplet which 
preserve the shape of the boundary ^ . 

Finally, to conclu de th is section, we reformulate the 
boundary condition (3.8) using the operator formalism 



in CFT. This is a purely technical step that will be used 
in the next section, when we analyze the consequences 
of (3.8) for the trial edge theory of the FQHE. It is a 



standard procedure in boundary CFT [301 SS] ■ The fields 
a{w) and (w) can be expanded in Fourier modes on the 



a{w) 



a' {w) 



E 



(3.10) 



^ e^"(a„)t (3.11) 



with the hermiticity condition (a_„)^ = (a^)„ (in this 
Euclidean field theory, the hermitian conjugate must be 
taken after continuation to real time on the cylinder, 
namely x — >■ it, so — —w). One has similar expan- 
sions for a{w) and a^{w). The boundary condition (3.8) 
can be written in terms of the modes as (a^)„ — a^n for 
any n € Z + ha- More precisely, this identity must hold 
while acting on a boundary state \B) 



\B) 



. 



(3.12) 



Such boundary states are known in the CFT literature 
as Ishibashi states JlOl H5] . It is convenient to think of 



the constraint (3.12) as the expression of an intertwiner 



between the chiral and the anti-chiral representations of 
A. Since the operator product expansions of a{w) and 
a^w) generate the full chiral algebra A, and because we 
are assuming that the representations of the chiral alge- 
bra A that appear in this paper are irreducible, Schur's 
lemma implies that the state \B) is completely fixed, up 
to a global normalization constant. 

Before we go ahead and analyze the consequences of 
these boundary CFT ideas, let us point out that other 
technical choices are possible for the analysis carried out 
here. We have adopted a grand-canonical formalism in 
order to be able to write the boundary condition (3.8 1, 
which violates neutrality. Neutrality is restored thanks 
to the fact that the particle number is not fixed. One 
could have adopted other conventions. For instance, one 
possibility would be to work in the canonical ensemble, 
and then focus on the neutral suhalgehra of the chiral al- 
gebra A, which is generated by all the neutral operators. 
For instance, in the U{1) sector, the neutral operators 
are the current idip{z) and its derivatives. In the statis- 
tics sector, the neutral subalgebra contains the stress- 
tensor, which modes generate a Virasoro algebra, and 
possibly some other local operators, which yield some 
extension of the Virasoro algebra (like a Wfe-algebra for 
Zfe-parafcrmions, see [33H3S])- Then one would have 



found a boundary condition analogous to (3.8) but for 



the neutral currents rather than for the operators a{w) 
and a(w) themselves. Also, another appealing possibility 
to circumvent the problems caused by the violation of 
charge neutrality, while working in the canonical ensem- 
ble, would be to use "shift operators" that would map 
the CFT vacuum with N charges, \N), onto the one with 
± 1 charges. These operators are not local. In terms 
of such a shift operator S, one would obtain a boundary 
condition of the form a{w) = S a^(u') S. The locality 
of the boundary condition would be somewhat hidden 
in this kind of expression. That is why, in this section. 
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rather than deahng with these shift operators, we decided 
to use a grand-canonical formulation in order to reach the 
local boundary condition (3.8), which should appear as 



more natural to the reader. Of course, although all these 
technical conventions need to be treated carefully for the 
global consistency of the argument, they will have no in- 
fluence on our final results. 



C. Back to the droplet 

In this section, we want to go back to the plane, where 
the particles fill a droplet of radius R. We want to un- 
derstand how one should handle the state 



1 



exp 



Jc 



when it appears inside correlators of the form 



|0) |0) (3.13) 



(7V|(7V|0(Ci)...(/.(Cp) 
1 



(3.14) 



A^ Z 



exp 



N 



|0)|0). 



Each of the operators (j){Cj) is one field a{(j), a{Q), 
'^HCj)i oi' ^HCj)- They are all lying outside the drop- 
plet: \Cj\ > R. We are interested in these correlators 
in the "scaling region" , which we now define, we first 
fix some number M > 0, and consider the correlators of 
the form (3.141 such that p < M. The non-zero con- 



tribution to the correlator (3.141 comes from the term 



generated by the exponential that has exactly the right 
total charge. This charge is contained in the interval 
[N - ^^,N + ^■]. Then we consider the limit iV -> oo, 
keeping M fixed. In that process, the radius R of the 
droplet grows, so one has to push the operators 4>{Cj) 
such that they stay out of the droplet (for example one 
can rescale the positions Q H> Q x R). In the scaling 
region, only terms with a number of particles within the 
range [N — + ^] matter. Different particle num- 

ber should in principle correspond to circular droplets 
with different radius R + SR, however we have defined 
the scaling region precisely such that SR/R — >■ when 
the number of particle goes to infinity, so the variations 
of the size of the droplet become negligeable. Therefore, 
in what follows the radius of the droplet is always R, even 
if the number inside it can fluctuate around N. 

Now we are ready to apply the ideas of the previ- 
ous section. If the screening hypothesis holds inside the 
droplet, then when the number of particles goes to infin- 
ity, the droplet becomes equivalent to a boundary con- 
dition at l^l = i? for the non-chiral CFT that remains 
outside the droplet. The exterior of the droplet can be 
mapped onto the right half-cylinder by the conformal 
mapping z i-^ w ~ ^log(z/it!), and we know that the 
boundary condition on the cylinder is (3.8). Again, the 
boundary condition (3.8) requires some fine-tuning of the 
parameter A, such that configurations with different par- 
ticle numbers contribute with equal weights. Let us skip 




R 



\B{R)) 



FIG. 4. In the thermodynamic limit, assuming screening, the 
state (3.131 is equivalent to the conformal boundary state 
\B{R)), in the so-called "scaling region". 



this detail for now. We reach the conclusion that, in ra- 
dial quantization in the plane, the boundary condition 
is again encoded by the boundary state (3.12) up to a 



scale transformation (in order for the boundary to be at 
l^l = i? rather than at \z\ = 1) 



\B{R)) = R' 



\B) 



(3.15) 



At this point, the reader might be worried by a techni- 
cal aspect: a neutralizing background was explicitly in- 
cluded in the previous section, while here we have traded 



the neutralizing background for the factor e 



V{z,z) 



the 



integration measure for the particles. However, since the 
boundary state \B) is already a sum over all the possible 
charge sectors with equal weights, this does not affect the 
expression of \B). In other words, when inserted in the 



correlators in the scaling region, the state (3.13) can be 
safely replaced by the boundary state \B{R)) in the limit 
N oo (up to a global normalization factor which still 
needs to be fixed). 



D. Consequence for the overlaps 
between the edge states 

Now we come back to the edge states which we defined 
in section [IIB[ and explore the consequences of screening 
for the overlaps between these wavef unctions. Following 
the formula (2.8), we define 



*(„|(zi, ..,zn) 



1 ^ 

— Y[aiz,)\0) (3.16) 



1=1 



where we have introduced the rescaling operator 1 . 
Her e, R is again the radius of the planar droplet defined 
by (3.4), and ALq = Lq — {N\ Lq \N) measures the con- 



formal dimension relative to the one of the vacuum {\N) 
is the vacuum with N charges). For an angular momen- 
tum eigenstate (Lq \v) = \v)), one can easily check 
that the factor l/R^^° ensures that ^f^^i (zi, . . . , zjv) is 
dimensionally homogeneous to the ground state wave- 
function ^(jv| (zi, . . . , zat), Zi and R being two lengths. 



The normalization (3.16) of the wavefunctions for the 



edge states will also allow us to express the result (3.23) 
in a particularly simple form. 



13 



Actually, the definition (3.16) is valid for the neutral 



edge excitations, and it is implicitly assumed that the out 
state {v\ is a state with charge N. To be able to express 
our final result ( 3.23 1 in a more general form, we also need 



to include charged excitations. Therefore, for a state {v\ 
with charge N + n {n can be positive or negative), we 
define the wavefunction for the excited state ^'^^| as 

(^1, ■ • • , Zn, Zn+1, ZN+n) (3.17) 




1=1 



where the coefficient A is the same as in (3.13). Now 



that our conventions for the edge states arc fixed, we 
can consider the overlap between two wavefunctions 
*(„j|(zi,...,2;Ar) and (zi, . . . , ^tv). 



(3.18) 



[^{vi\{zi, ■ ■ ■,Zn)]* • ■ -iZn)- 

This is an overlap between two neutral edge excitations. 
More generally, between two charged excitations with 
some charge N + n, the overlap is defined as 



1 



N+n 



{N + ny. 

[^{vi\{zi, . ■ . , ZN+n)] * *(i,2| (^1' 



•(fz, (3.19) 

. . . , Zjv+n)- 



The overlap between two wavefunctions with different 
particle numbers is always zero. Using our definition 
( 3.16||3.17 l, these overlaps are equal to 



1 



{V1\{V2\R 



-ALn-ALn 



A / e^(^'^)d2za(z)a(z) 



(3.20) 

\o)W)- 



According to the previous section, the whole expression 
in the second line can be replaced with the boundary 
state R^o+Lo 1^^^ least as far as we are in the scaling 
region. This means that the charge of the states {vi\ 
and {v2\ must be kept to some value N + n, where \n\ < 
M, and M is fixed while we send N to infinity. This is 
precisely what we do here. Then, we have 



JV-s-oo 



z 



N 



v,\{v2\ \B) (3.21) 



with /lAT + ZiAT = {N\Lq\N) + {N\Lo\N). Using some ba- 
sis of states an^an^ . . . (a^ )mi (Q ^)m^ • • ■ |0), one can easily 
show that the property ( |3.12[ ) implies that for any \vi) 
and 1 1^2) 



{vi\{v2\\B) — (w2 X constant 



(3.22) 



where the constant does not depend on \vi) and \v2), 
and comes only from the global undetermination of the 



normalization of \B). The property (3.22) is 
known property of Ishibashi states (see [40^ 
ally, it could even be used as a definition of 



state, instead of (3.12) 



The normalization of \B) is fixed by the 
that the ground state wavefunction 5'(Ar| is 



a very well- 
|43]). Actu- 
an Ishibashi 

requirement 
normalized: 



be equal to Zn/R''" 



Thus the constant in (|3.22|) must 
The final result is 



then 



N- 



{V2 \vi) 



(3.23) 



This formula, which is a direct consequence of the gen- 
eralized screening hypothesis, is the main result of this 
paper. It shows that the linear mapping (3.16) from the 



(dual of the) Hilbert space of the CFT to the space of the 
physical edge states becomes an isometric isomorphism 
in the thermodynamic limit. This is a precise formula- 
tion of the long expected "bulk/edge correspondence". 
It implies, in particular, that the underlying CFT has a 
positive definite inner product, or in other words, that it 
is a unitary CFT. In the MR construction (section IC), 
only the rationality of the CFT for the statistics sector 
was assumed. We see th at, if generalized screening holds, 
then we arrive at (3.23), which implies that the CFT is 
unitary. This is one additional argument that shows that 
the use of a non-unitary CFT in the statistics sector can- 
not be consistent with the fact that all the connected cor- 
relations of local operators in the bulk are short-range, 
as it clearly leads to contradictions (for previous argu- 
ments, see Refs. [Ml Hi])- Therefore, for consistency, the 
FQHE trial states given by the MR construction should 
correspond to rational and unitary CFTs only (for the 
statistics sector). 



The formula (3.23) will also play a key role when we 



study the entanglement spectrum in part |IV[ We will 
provide some direct numerical checks of this result in 
section ITlIFI In the next section we derive another result 
that is directly related to this bulk/edge correspondence. 

To conclude this section, let us comment on the factor 
A that appears each time one has to deal with particle 
numbers that differ from N. We have for instance 



(JV+ll |*(Ar+l| 



1)) 



N- 



1, 



(3.24) 



while by using the definitions (3.2), (3.17) and (3.19) 



(Af+l| |*(Ar+l| 



1)) 



A 



Z 



N+l 



/^'iN + l+'lN + l Zj'^ 

(3.25) 

In the previous sections, we explained that the coefficient 
A had to be tuned such that it gives rise to the result 
(3.23) for charged excitations (not only for the wavefunc- 
tions for neutral excitations). Thus, the formula (3.24) 



is rather a definition of the coefficient A than an actual 
result. Indeed, in general A may depend on the radius 
R, and therefore on the number of particles N. How- 
ever, once A is fixed, there is no other free parameter in 

{N+n\ |*(Ar+n|)) , 



(3.23). For instance, by evaluating ((5' 
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one finds that Zt^+nl^N , 
where the coefficient A is no longer a free parameter. As 
an exercise, one can check that this is consistent with di- 
rect free-fermion calculations in the integer quantum Hall 
effect, or with the results of the semi-classical expansion 
of Wiegmann and Zabrodin for the Laughlin wavefunc- 
tion EHl. 



E. Equality of edge and bulk CFT correlators 

In this section we show that screening, reformulated 
as the conformal boundary condition (3.8), implies that 



the equal-time correlators measured along the edge of 
the quantum Hall system are equal to the correlators 
in the CFT that is used to construct the ground-state 
wavefunction (1.6). For instance, we can compute the 
correlator of particle creation/annihilation operators cj S 
along the boundary of the droplet (|zj| = \z'j\ — R), 

{{^\cHz[)...cHz'Jc{z,)...ciz,,)m. (3.26) 



The ground state | 4')) is the wavefunction (1.6) with N 
particles. In first quantization, the correlator (3.261 is 



A 



(N-nV. Jc 

[5'(zJ,...,Z^j,Z„+i. 



(3.27) 



, ZAr)]**(zi,...,2Ar) 



A (iV|a(2'i) . . . a{z'^) (7V| a(zi) . . . a(z„) x 
1 



■ exp 



A y e^(^^^)rf2^ a{z)a{z) 



|0) |0) 



A A" {N\a{z\) . . . a(zO (A^ | a{z^) . . . a(z„) \B{R)) 



The factor A is the product f]^^^ e^^^^^^f )/2e^(S'4)/2. 
Since we assume that Viz^ z) is rotationally invariant 
and the Zp and z^ are along the circle of radius R, one 
has A — e"^(^). In the last line we have replaced 



the state (3.13) by the boundary state (3.15), as ex- 
plained in section [ill C| Finally, we use the fact that the 
boundary state \B{R)) implements the conformal bound- 
ary condition a^{w) — a{w) on the cylinder, which can 
be mapped onto the plane by the conformal mapping 
K7 > z = i? exp . This leads to the boundary condi- 
tion along the circle of radius R in the plane 



(3.28) 



The factors , J§ appear because the operators a and 
a transform covariantly under conformal transformations 
(see [27H29] ). These two factors are equal to ^ and 2|2 

i?is 



respectively. Thus, the boundary condition at |z| 

(z/z)^" 



z) (recall that ha = ha). In the end. 



the correlator ( |3.26 ) converges to the following correlator 



in the chiral CFT in the plane 



n 



p=i 



X {a^{z[) . . . a^(4)a(zi) . . . a{zp)) . (3.29) 

In particular, the particle propagator along the edge is 
proportional to the two-point correlator in the chiral 
CFT, {a^z')a{z)) 



{{^\cHz')ciz)m (X — 



(3.30) 



in the thermodynamic limit. This shows that equal-time 
correlators evaluated along the edge are given by cor- 
relators in the CFT that is used to construct the trial 
wavefunction for the ground state (1.6). This result has 



been assumed in many places in the literature, although 
no general argument has ever been given. It generalizes 
the one obtained by Wen for the Laughlin wavefunction 
in Ref. [T5]. For recent work on this topic, see also [15] . 
Note that we have restricted the result to equal-time cor- 
relators because we do not address physical Hamiltonians 
in this paper (see, however, [SUl \5T\ ). 

The discussion in this section can be directly extended 
to the case of equal-time correlation functions of quasi- 
particle and quasi-hole operators along the edge. In the 
thermodynamic and assuming screening, the same calcu- 
lation base on the boundary condition (3.8) can be done. 



leading to the equality (up to normalization and phase 
factors) between these correlation functions and the cor- 
relators in the bulk chiral CFT that underlies the trial 
wavefunction. 



Numerical checks 



The formula (3.23) can be checked numerically. In this 



section we present numerical evidence that shows that it 
holds for the Laughlin wavefunction and for the Moore- 
Read (Pfaffian) wavefunction. In both cases, we do a 
Monte-Carlo (MC) simulation for a system of N particles, 
which is tractable both for Laughlin and Pfaffian states. 
The MC allows us to estimate numerically the ground 
state expectation value of any observable 0{{zi}) that 
depends on the positions z^'s 



{{^\o{{z,})m. 



(3.31) 



1. Laughlin at u — 1/3 

The edge states for the Laughlin wavefunction are 
given explicitly in section [IIB| We have 



{N\Jk^...Jkp 



(3.32) 
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where 



2. Moore-Read (Pfaffian) at u = 1/2 



Sk 



N 

E 



(3.33) 



and SI is the complex conjugate of Sk- It is the right- 



hand side of ( 3.32 1 that we measure numericahy with MC 



techniques. The result predicted by (3.231 when N ^ oo 
is 

{{^\Sl...SlSk'^...Sk'J^)) 

[jk'^...Jk'^J^k,...J-k,). (3.34) 
For jy = 1/3, iV = 100, with 10^ MC steps and with the 



quadratic potential (1.3) corresponding to the plane, we 
find the following results for the first exited states. They 
are in very good agreement with our analytic prediction 
in the N ^ oo limit. 





MC 


analytic 




1.007 


(JiJ-i) = 1 


{{^\s*^s2m 


2.017 


(J2J-2) = 2 


{{^lisirs^m 


0.003 


(j2Jii) = 


{{^\s*sisiS2m 


2.034 


(J2J1J-1J-2) - 2 


{{^\{s*,)HS2?m 


8.048 





For the Pfaffian state, there are two types of edge exci- 
tations: the excitations in the C/(l) charge sector, which 
are similar to the ones of the Laughlin state, and the 
excitations in the Majorana sector. For the U{1) excita- 
tions, we find for v = 1/2, N = 100 and with 10*^ MC 
steps: 





MC 


analytic 


{{^\sts,m 


1.002 


(JiJ-i) = 1 


{{^isis^m 


1.993 


(J2J-2) = 2 


{{^\{st)'S2m 


0.005 


(JzJ^i) = 


{{^\{st)Hsirm 


1.995 


(Jfj'-i) = 2 


{{^\{s*2?{S2?m 


7.934 


{JlJ'-2) = 8 



In the Majorana sector, the excitations are of the form 
(2.14). For instance, with two excited fermion modes 



only, we get the overlaps 



(3.35) 



where i^ni,n2 the following ratio: 



Pf 



n2~ i 



-1 
Z1-Z2 



"N 



Rni+n2 X Pf 



Z1-Z2 





<i,j<N 



-N 



ZN-1~Zn 





(3.36) 



Similar formulas hold for more fermionic excitations, for example for i^ni,n2,n3,n4j etc. Again, the right-hand side of 
(3.35) can be measured numerically in a MC simulation. When N ^ 00, we expect 



(3.37) 



We have checked this for a few matrix elements for sizes N = 10, 20, 30 and 100 (each of them with 10^ MC steps). 
The results are in good agreement with our analytic prediction, although the convergence is slower than in the U{1) 
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sector. 





TV = 10 


iV = 20 


iV = 30 


N = 100 


analytic {N oo) 




1.371 


1.288 


1.231 


1.124 




2 ' 2 2 ' 2 


1.460 


1.323 


1.268 


1.151 


(ihi'disib 5ih i) = l 

\ 2 2 2 2/ 


{{^\Ft 5^3,1 1*)) 
2 ' 2 2 > 2 


1.432 


1.393 


1.331 


1.189 


( lb 3 lb lb 5 lb 3) = 1 

\ 2 2 2 2/ 


2 ' 2 2 ' 2 


-0.033 


-0.004 


-0.003 


-0.001 


(^|^5^„5V'_l) = 


{{^^J\Ft 3 5 zFi.3.|.||*)) 

2'2'2'2 2'2'2'2 


1.617 


1.764 


1.700 


1.381 


( tb 1 lb 3 lb 5 itmb iib sih sib i ) = 1 

\ 2 2 2 2 2 2 2 2/ 



IV. ENTANGLEMENT SPECTRUM 

We are ready to apply the formalism developed in parts 
|TT] and |III| to the analysis of the entanglement spectrum 
of the ground state wavefunction. In particular, the uni- 
versal result (3.23); which relates the overlaps between 



the edge states to the inner product in the CFT space, 
will play a key role in our analysis. Our main focus is on 



real-space partition (RSP, see section ID I. However, the 
techniques we use can also be applied to the case of the 
particle partition (PP) in t he so- called "scaling region", 
as explained below (section IV D I . 



In this part it is more convenient to use a second quan- 
tized formalism. By definition, in second quantization. 



the wavefunction (1.6) becomes 



(4.1) 



^izi,...,ZN)c\zi)...c\zN)m 

where the c^(z)'s create fermions/bosons at position 
z and 1 0)) is the fermionic/bosonic vacuum annihi- 
lated by all the c(z)'s. The modes obey the canonical 
(anti-)commutation relations {c{z), c^(z')} = S^^''{z — z') 
(fermions) or [c{z),c{z')] = S^'^\z — z') (bosons). To de- 
fine the RSP, we fix i? > 0, and we define A as the disc 
of radius i?, A = {z G C; \z\ < R}, and B as the comple- 
mentary subset B ^ C \ A. Our goal is to compute the 
Schmidt decomposition 



N 

EE' 

Na=0 k 



-UNA,k)/2 



where {| (as well as {| ^^^^k))}) is a set of 

orthogonal states with Na particles in part A (and Nb = 
N — Na particles in part B). The set of pseudo-energies 
£,{NA,k) is the entanglement spectrum [131. 

Since the bipartition AUB is rotationally invariant, the 
angular momentum in part A, Lf, is a good quantum 
number. The Schmidt eigenvalues/eigenvectors can be 



ID 



classified according to Lf. As claimed in section 
general, there is a non-degenerate lowest pseudo-energy 
^ at some values Nao and L^q. These values depend on 



(ze) J2\NA,k){NA,k\R-^<' 




FIG. 5. We use the completeness relation (4.41 to decompose 
the conformal correlator (]~[^ a{zi)^ as a sum of products of 
one correlator which involve only particles in part ^4 (1^1 < R) 
with another one which involves particles in part B {\z\ > R). 



the radius R. In general, the radiu s R is the one of a 
droplet with Nao particles given by (3.4) 



and is the 



angular momentum of the ground state with A'^o 

particles. We define A^, AA^^i and AL^ by subtracting 
of these values. 



A. Decomposition of the ground state from 
a completeness relation in the CFT space 

In this section, we use a completeness relation to de- 
compose the ground state as a sum of products of terms 
of the form | ^P^)) | ^^))- In the end, these terms will 
turn out to be orthogonal only in the thermodynamic 
limit, when Na,Nb — >■ oo (and assuming screening), so 
this will be a true Schmidt decomposition only in that 
limit. However, it will not be difficult, once our proce- 
dure has been exposed, to adapt it to take into account 
corrections to scaling. This will be explained later in sec- 
tion |IV^ For now, we initiate the process of calculating 
the decomposition in a naive way, relying only on the fact 
that the ground state wavefunction is a conformal block. 

Like in section |IIB[ we think of the CFT in radial 
quantization, and we use a complete set of states in the 
CFT Hilbert space to write the identity operator as 



\NA,k) {NA,k\ 

NA,k 



(4.3) 
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In this expression, Na is the number of particles, and 
k stands for all the excited states in a given Na sector 
(called the descendants in the CFT literature). Strictly 
speaking, in radial quantization, the operator (4.3) ac- 



tually acts on the CFT states corresponding to the field 
configurations on a circle at radius \z\ — 1, so we need 
to rescale this operator such that it acts at \z\ = R (see 
also Fig. m 

1 = ^ \Na, k) {Na, k\ i?-^" . (4.4) 

NA,k 



Now we write the state (4.1) as 



At vz 



(4.5) 



N 



(iV|7^exp 



mm 



Note that |0) is the CFT vacuum, while |0)) is the 
physical fermionic/bosonic vacuum (annihilated by the 
c(z)'s). The symbol TZ denotes radial ordering (see [27l - 
[5U]). When one expands the exponential, only the term 
with N particles remains because of the projection onto 
the vacuum with N charges |iV). This expression is then 
exactly equal to (4.1 1. The integration over C = AU B 



can be split into an integration of A and another over B, 
so the exponential is equal to 



7?,exp 



IB 

X 7?.exp 



A / e^'-^'^^/^d^z a{z)^c\z) 



A / e^^^'^^/^d^^ a{z)(g)c^{z) 

. J A 



Then we insert the completeness relation (4.4 1 between 
these two exponentials (Fig. [5]). The result can be writ- 
ten as 



Z 



N 



Na o '^Nbo 



JVa=0 k 



A 



^\NA,k) 



where Nbo = N — Nao, and 



{NaM 

72.exp 

^\NA,k) 

TZexp 



A- 



{NA,k\ 



RALo 



A / e^^^'^^/^d^z a{z)(E)cHz) 

J A 

^ {N\ X 



|0) |0)) 



(4.7) 

(4.8a) 
(4.8b) 



A / e^("'")/2^2^ a{z)®c^z) 

J B 



R'^^°\NA,k) |0)) 



We have used the notation ALq = Lq — {Nao \ Lq \Nao)- 
The normalizing factors Zj^^^^ and Z^siga are fixed such 



that the states 



^{Nao\ 



and 



correspond- 



ing to the vacuum |iV4o) are both normalized 



-, „ Nao Nao 

- ^ / n ^""^'-^^d^^, I {nao\ n io) 

AO- J A 



Znao - 



(4.9a) 



Nb 



^0- "'s i=l i=l 

(4.9b) 



The states 



^{NaM 



appearing in the decomposition 
(|4.7|) should be viewed as the edge excited states that we 
II 



constructed in part 



^\NA,k) 



Similarly, the states 

can be viewed as the edge excited states for part B (al- 
though, strictly speaking, we have not constructed these 
states in section |IIB[ the extension of our formalism to 
the part B is straightforward). One difference with the 
edge states introduced previously is the integration do- 
main: here the positions of the particles are restricted to 

the domain A, so the 4*^^ '^l))'^ ^'"^ projections onto 
the domain A of the edge states of part |TT] With that 
minor distinction in mind, we see that the formula (4.7) 



is a natural decomposition of the ground state as a sum 
of products of edge excited states for part A and part 
B. It is quite close to the Schmidt decomposition we are 
looking for, except that the states appearing in the right- 
hand side of (4.7) are not necessarily orthogonal. In the 



next section we argue that, in the "scaling region" (to be 
defined below), the states are othogonal if the screening 
assumption holds, such that we can use our result for the 



(4.6) overlaps between the edge states (3.23) 



It is worth emphasizing that, although the decomposi- 
tion (4.7) may look quite formal, it becomes very explicit 



when one considers concrete examples. In the case of the 
Laughlin wave function, for instance, the decomposition 
of the ground state in terms of the edge state wavcfunc- 
tions becomes very elegant; it can be found explicitly in 
[TS]. For the Moore-Read (Pfaffian) wavefunction (2.12), 



it is an instructive exercise to write down the decompo- 
sition of the Pfaffian in terms of the edge excited states 
(|2T4|. 



B. 



Pseudo-energies at the leading order 
in the scaling region 



Let us start by defining what we call the scaling re- 
gion. It is the extension of the concept of "scaling re- 
gion" that we used previously in part |IIH which needs to 
be adapted to the bipartition AU B. Roughly speaking, 
it is the set of Schmidt eigenvalues and eigenvectors for 
which ANa is small compared to N, and ALf is small 
compared to Lf^. More precisely, we imagine that the 
radius R of the disc A depends on N such that Nao/N is 
fixed. This is the situation that is usually considered in 
numerical studies, where A and B correspond to the two 
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hemispheres of the sphere after stereographic projection, 
with a ratio Nao/N = 1/2 (for an even number of par- 
ticles on the sphere). Then we fix some positive number 
M, and keep all the eigenvalues and eigenvectors such 
that \ANA\,\ALf\ < M. Finally, we take the thermo- 
dynamic limit N oo. This scaling region is the one 
where the entanglement spectrum is conjectured to pos- 
sess the scaling property [TS]. Note that, by definition, 
in the scaling region, both Na and Nb go to infinity. 
In this section we analyze the limiting behavior of the 
pseudo-energies in this range of quantum numbers AA^^, 
ALf. We do this only at the leading order; subleading 
corrections will be tackled in the next section. 

First, we focus on the prefactor appearing in the right- 
hand side of (4.7 1. This is a ratio of partition functions 



result (3.231 



Z 



Na o ^Nbo 



e^V[-UNAo+fNso-lN)l'A (4.10) 



where Jnao ~ ~ ^^g ^at^q , etc. Each of these three terms 
in the exponential is the free energy of a statistical sys- 
tem of itinerant particles in part part _B, or in C. To 
discuss them, it is natural to adopt a statistical mechan- 
ics point of view as in part |III[ The following argument 
is standard in that context; it has been used in several 
works (see for instance ^52^) to support the area law for 
the entanglement entropy |12| or other similar quanti- 
ties such as the (logarithm of the) ratio (4.10) [53 . The 
leading contribution to the free energy is a bulk part, 
which in general could be expressed as a functional of 
the local density of particles inside the domain. This 
bulk part is usually extensive for statistical systems with 
short-range interactions, or super-extensive if the inter- 
actions are long-range, like Coulomb interactions in the 
plasma mapping for the Laughlin wavefunction. Another 
contribution, which is subleading, comes from the bound- 
ary of the statistical system. It usually scales linearly 
with the boundary, and may be interpreted as a surface 
tension. In particular, the free energies /jv^q and fNeo 
should both contain a surface tension proportional to the 
length of the boundary between A and B (that is 2ttR). 
The bulk parts of Jnao ^'^'i Inbo cancel the bulk part of 
/at , leaving only the surface tension term at the boundary 
between A and B. Therefore, we are left with 



log 



^Nao^Nbo 

Zn 



a 27ri?-7-fo(l) (4.11) 



when R is large compared to mean particle spacing close 
to the cut (which plays the role of the UV cutoff). The 
coefficient a is non-universal and depends on the UV cut- 
off. We have also included a possible 0{1) term, noted 
—7. This term is dimensionless, and it is universal: it 
is the topological entanglement entropy jlOl [TT] . The 
presence of this term here can be shown by adapting the 
original argument of KP; this is explained in more details 
in the appendix \K\ 

Second, assuming screening, in the limit Na — > 00, the 



overlaps between the states 



^{NaM 



are given by our 



^fNA 



Na 



{NA,k'\NA,k} = Sk,k' 



(4.12) 

Here we have used the fact that the basis of states 
{\NA,k)} in (4.31 is orthonormal. States with different 
particle number N'j^ ^ Na are always orthogonal. The 
same formula can be obtained for the overlaps between 
the states in part _B by a straightforward extension of the 
arguments developed in part |III| 



^\NA,k) 



^\NA.k') 



Nb 



{NA,k\NA,k'} = Sk,k' 



(4.13) 

One might be worried about the integration domains: 
here the positions of the particles Zi are integrated only 



over part A (or B), while the result ( 3.23 ) was claimed for 
overlaps computed by integrating over the whole plane C. 
This does not make any difference though, as long as we 



focus on the states 



^\NA,k) 



or 



^{NaM 



that are in 



the scaling region. For these states, the contribution of 
the integration over the whole plane is exponentially sup- 
pressed by the factor e^'^'^^ outside the droplet of radius 
R, so the difference between the overlaps computed by in- 
tegrating over A (B) or over C is exponentially small. The 
choice of the non-universal coefficient A in (4.81, which 



must be tuned in order to give the correct weight to the 
terms with different particle numbers, might be sensi- 
tive to fact that one is integrating over A 01 B rather 
than the whole complex plane C. However, as discussed 
at the end of sction |IIID[ there is always one choice of 
A that leads to the conformal boundary condition (3.8) 



and to the universal result (3.23). Therefore the limits 



(4.12) and (4.13) are correct, and the subtelties about 



the integration domains will become relevant only in the 
next section, when we deal with subleading corrections 
to these formulas. 



Let us come back to the decomposition (4.7). We have 



reached the conclusion that, in the scaling region, the 
two sets of vectors in the right-hand side are orthogonal 
at th leading order. Thus, in that limit, (4.7) is really a 
Schmidt decomposition. In particular, this implies that 
in each Lf subsector, the Schmidt rank matches exactly 
the dimension of the corresponding subspace of angu- 
lar momentum eigenstates in the conformal field theory. 
Also, the Schmidt eigenvectors turn out to be exactly the 
edge states constructed in section [Hj up to the subtelties 
involving the integration domain. The pseudo-energies 
are all equal to 



^{ANa, k) ^ a 2ttR - 7 + o(l) 



(4.14) 



in the scaling region, independantly on ANa and on 
the descendant labelled by k. The splitting between the 
pseudo-energies appears only at higher order, when one 
tackles subleading corrections to the formulas ( 4.12|4.13 1. 
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C. Corrections to scaling, RG analysis, and locality 
of the pseudo-Hamiltonian 

In this section we use ideas that come from the field of 
surface critical phenomena (for a review, see Ref. |54| ) . 
Like bulk critical phenomena, they can be understood 
with the renormalization group (RG). Let us consider a 
system which is critical in the bulk, such as a critical Ising 
model in d dimensions. At the surface (which is d — 1 di- 
mensional), let us imagine that the spins are free. Let us 
imagine that one can turn on a magnetic field at the sur- 
face. Then the spins at the surface tend to align with the 
magnetic field. Thus, the surface of the system under- 
goes a transition, although the bulk is still critical. Under 
the RG flow, the surface of the system flows towards a 
fixed where all the spins are aligned. This is called a 
boundary RG flow; which brings the system from one 
type of boundary condition to another one. A boundary 
RG fixed point is a scale-invariant boundary condition. 
For most systems, scale-invariance extends to conformal 
invariance, and we get a conformal boundary condition 
(which means a conformally invariant boundary condi- 
tion) . At the vicinity of a conformal boundary condition, 
the scaling behaviour can be understood in terms of per- 
turbing operators along the boundary of the system. In 
our example of the Ising model, the operator which is 
coupled to the magnetic field at the surface is the one 
corresponding to the local magnetization. This operator, 
which we note (f>a{x), is located at the surface, therefore 
it is called a boundary operator. Like in the more famil- 
iar case of bulk critical phenomena, one can classify the 
boundary operators as relevant, irrelevant, or marginal, 
depending on their scaling dimension ha- When such an 
operator appears as a perturbation at the boundary, it 
adds a term of the form Aa / (i'^^^X(/>a(x) to the action of 
the theory. The coupling Aa scales like ^''i+i-'*^ where 
£ is some UV cutoff, such as the lattice spacing if our 
system is a statistical model on a lattice. Under the RG 
flow, 0a is said to be relevant when ha < d—1, irrelevant 
ii ha > d—1, and marginal ii ha = d—l. Generically, all 
the operators which respect the symmetries of the sys- 
tem are expected to appear as a boundary perturbation. 
For more information on surface critical phenomena and 
boundary RG flows, we refer the reader to [MlESj. Now, 
let us explain how to use this framework for analyzing 
the ES of wavefunctions given by conformal blocks. 



1. Locality 



Our main result (3.23) and its variations ( 4.12||4.13 ) 
that are useful for the calculation of the entanglement 
spectrum, rely on the fact that, in the thermodynamic 
limit and assuming generalized screening, we are left with 
a non-chiral CFT that lives outside the droplet, con- 
strained by a boundary condition along the edge of the 
droplet. The interior of the droplet decouples from the 
exterior thanks to screening. In this framework, it is nat- 



ural to include boundary perturbations that modify the 
conformal boundary condition. The action of the field 
theory, 5cft, is then modified along the circle |z| = i? 
by boundary perturbations S'cft — > S'cft + St- The 
latter are of the form 



Sb — ^ Aa 



\dz\ 4>a{z) 



(4.15) 



\z\=R 



where the (jja^s are some local boundary operators with 
scaling dimensions /la, and with coupling constants Aa. 
The boundary condition which we discussed in part III 
should not be stable under the RG flow, which means 
that there can be no relevant perturbation, namely all 
the scaling dimensions satisfy /la > 1- The boundary per- 



turbation (4.151 might look completely generic, however 



one should emphasize that it implicitly assumes locality, 
in the sense that it is a sum of local operators. This holds 
thanks to the locality of the action of the CFT outside the 
droplet, and thanks to the screening assumption inside 
it. Without screening, the degrees of freedom along the 
edge might be coupled at long distances through the bulk, 
which would typically lead to non-local perturbations of 
the action S'cft- Since we are assu ming s creening, this 
cannot happen, and the perturbation (4.15) has the most 



generic form. The coupling Aa scales with some power of 
the UV cutoff, which is of the order of the mean particle 
spacing close to the edge, or equivalently y^^^ where 
po is the mean particle density close to the edge. Thus, 



The perturbation (4.15) modifles our formula for the 



overlaps between the edge states in the scaling region 
(which implies y^~^/i? — >■ 0) 



(■"2 I 



{V2\e 



"^^ \vi) 



(4.16) 



Of course, the leading order in this formula is noth- 
ing but the universal result (3.231, however this new 



expression generates corrections to scaling. The flrst 
correction comes from the least irrelevant operator 
{i.e. with the smaller ha) and leads to a term of 
order (^p^" V^)''"^^ ~ {l/y/NAo)''''-\ The de- 



nominator {e-'"^") = {NaoIb-^" \Naq) is fixed by the 



is normalized: 



requirement that the state 

((*^.o| l^fiv.ol)) = 1- By redefining Sb ^ Sb + 
constant, one can absorb this denominator in the defi- 
nition of Sb itself. This is what we do in the following, 
which yields 



("2 1 



{V2\e 



(4.17) 



and the same formula holds for part B. However, since 
the boundary perturbation Sb is not necessarily the same 
for A and B, in what follows will label it by Sb{A) or 
Sb{B). 
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Now let us come back to the problem of computing the 
entanglement spectrum. The operator 



(4.18) 



acts in the Hilbert space of the chiral CFT. Note 
that, because Sb{A) and Sb{B) are both given by inte- 
grals of local operators, the commutators [Sb{A), Sb{B)], 
[Sb{A), [Sb{A), Sb{B)]], etc. are also given by integrals of 
local operators. This can be easily shown by writing the 
products of several fields <i>a^ (zi), (22), . . . in ( |4.15 ) 



with the help of the operator product expansion. Then, 
using the Baker-Campbell-Hausdorff formula, we see that 
the operator 5'es in (4.18) is itself the integral of some 



local operators, of the form (4.15) 



The locality of the operator Ses is another important 
result of this paper. Like our first important result, which 
was the isometric isomorphism (3.23) between the phys- 



ical Hilbert space of edge states and the CFT space, 
the locality of ^es is a consequence of the generalized 
screening assumption, for wavefunctions in the MR con- 
struction. The reason why the locality of S'es is impor- 
tant is that this operator is nothing but the "pseudo- 
Hamiltonian" we are looking for. Its spectrum is the 
ES, up to the global additive constant (4.14). The cal- 



culation of the ES in the scaling region has thus boiled 
down to a standard renormalization group (RG) discus- 
sion: one needs to determine which local operators (j)a{z) 
(which can be relevant, irrelevant or marginal, depending 
on their scaling dimension ha) do or do not contribute to 
the expression of S'es (of the form (4.15 )) relying on sym- 
metry arguments. We will come back to this below; for 
now, let us explain why the ES is related to the operator 
S'es- 

We consider the singular value decomposition (SVD) 
or e 2 



AS(W^,fc) 



NA^k 



\UNA,k) [VNaM 



(4.19) 



By definition of the SVD, {\uNA,k)} and {\uNA,k)} are 
two sets of orthogonal vectors in the CFT Hilbert space. 
Here the A^(iV^,fc)'s are simply the eigenvalues of S'es- 
Next, we write the identity operator as 



1 = 



(4.20) 



\NA,k J 

and we insert this in (|4.6|), like we did previously in sec- 
rhis timi 

Znao ^Nbo 



tion IV A This time, however, we get 

A;(iv^,fc) 



Na k 



e 2 ^uMj, 



(4.21) 



The difference between this last identity and (4.7) is that 



the states in the right-hand side are now orthogonal, be- 
cause of (4.171. The operators Sb{A) and Sb{B) are also 



assumed to be self-adjoint here; this reflects the fact that 
they represent boundary perturbations of the action of 
the field theory, which must be real. Then, (4.21 ) is truly 
a Schmidt decomposition. The pseudo-energies are di- 
rectly given by the eigenvalues of S es (no ted A^{Na, k)), 
up to the global additive constant (4.14) 



^Na, fc) = a 2^i? - 7 + A^Na, k) . 



(4.22) 



This completes our proof of the "scaling property" con- 
jectured in il5i , which states that the entanglement spec- 
trum in the scaling region is the spectrum of a sum of 
local operators along the cut between A and B. Our line 
of arguments relies heavily on the generalized screening 
hypothesis, and makes use of standard RG arguments 
and field-theoretic tools. 



2. Pseudo-Hamiltonian, Hamiltonian of the 1 -I- Id CFT, 
and role of the stress-tensor 

So far, we have shown that when the generalized 
screening assumption holds, the pscudo-Hamiltonian S'es 
is the integral of a local operator in the 1 + Id CFT 
that underlies the wavefunction. We have not deter- 
mined what this operator is yet. We are now going to 
argue that, in generic situations, and for real-space parti- 
tion (RSP), this operator is the Hamiltonian of the CFT, 
namely j^Lo for some "velocity" v, with additional per- 
turbations. Thus, the real-space ES of the trial wave- 
functions is the spectrum of a perturbed CFT. 

Once again, we look at the non-chiral CFT that 
lives outside of the droplet, with a perturbed conformal 
boundary condition along the boundary. In the case of 
our FQHE states, there is always a marginal operator, 
which is U(l) current idip{z), with scaling dimension 1. 
However, as a perturbing boundary operator, one gets 
the integral of the current along the cut, which is noth- 
ing but the number operator, so it plays a role only when 
sectors with different particle numbers are involved. The 
weights of the different sectors may eventually be fixed 
such that the spectrum does not depend directly on the 
number operator. Actually, this is exactly what we did in 
part |III| when we explained that the coefficient A could be 
chosen arbitrarily, and we tuned it such that the isometry 
(3.23) holds not only for neutral excitations, but also for 
charged excitations. We see now that the fact that this 
coefficient required some fine-tuning was simply due to 
te presence of a marginal operator. That being said, the 
role played by the U{1) current is very easy to take into 
account, so we can safely turn to the next boundary per- 
turbations, which are usually irrelevant. In particular, 
the stress-tensor T{z) appears generically as a perturba- 
tion. The coupling At has the dimension of a length, and 
is of order At ^ \/Po~^ ^ ^b- In statistical mechanics, 
the length 27rAT is known as the extrapolation length and 
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is ubiquitous in the study of surface critical phenomena 
|54j . Ahhough this has the dimension of a length, we 
use the notation v — AttXx, because v is going to be the 
velocity that appears in the ES. Generically, there are of 
course other perturbations, but for the sake of simplicity, 
let us assume first that this one is the least irrelevant one 
(this is correct for many trial states, but in principle, this 
assumption needs to be discussed case by case). Then, 
by restricting the analysis to this perturbation only, one 
finds that (3.231 is modified as follows 



(vi e 



«2 



(4.23) 



We assume that this formula is valid both for part A 
and B. Then, from the results of the previous section, 
we see that the entanglement spectrum is given by the 
eigenvalues of the operator 



S 



ES 



SbiA) 



(4.24) 



This means that the entanglement spectrum is the spec- 
trum of ^Lq in the chiral CFT that underlies the trial 
wavefunction, with a linear dispersion relation. Our re- 
sults thus support the claim made in [101 ES] . 

More generally, however, other boundary perturba- 
tions are present in (4.15 1, and may consequently appear 



in the ES. The latter is thus given by the spectrum of ^Lq 
in a perturbed CFT. All the perturbations that preserve 
the symmetries should generically appear, and this needs 
to be discussed case by case. For most trial states, we ex- 
pect the other perturbing operators ^a{z) appearing in 
S'es to be more irrelevant than the stress-tensor, namely 
they should have scaling dimensions ha > = 2. In 
that case the ES in the scaling region collapses onto a 
branch of excitations (a straight line) with velocity v as 
iV — oo. The corrections to scaling due to the more ir- 
relevant operators introduce some splitting between the 
pseudo-energies. We would like to warn the reader about 
the fact that, although the stress-tensor is itself irrelevant 
in the RG sense {hr = 2 > 1), it is crucial to take it into 
account to understand the structure of the ES, at least 
beyond the leading order where all the pseudo-energies 
are degenerate. When discussing the influence of the dif- 
ferent perturbing operators in S'es: most of them turn 
out to be irrelevant in the sense that ha > 1, but they 
can still have spectacular effects on the ES if they are 
less or equally irrelevant to T{z) (that is, /lo < 2). The 
standard RG terminology (relevant /irrelevant/marginal) 
must thus be used cautiously here, as we refer to the 
relevance of the boundary perturbations at the boundary 
of a 2c? theory, which differs from the relevance of bulk 
perturbations of the 1 -|- Id theory that gives the ES. 

It is also important to point out that, for most states 
in the MR construction, there are perturbating operators 
in the ?7(l)-charge sector, in the statistics sector, and 
possibly also mixed terms between the charge and statis- 
tics sectors. In particular, the stress-tensor of the full 
theory breaks down into two pieces: one for the charge 



sector r^(i)(2) and another one for the statistics sector 
T^{z). These two operators do not necessarily appear 
with the same coefficient Ay. Thus, when the statistics 
sector is non-trivial, like for the MR (Pfaffian) wavefunc- 
tion or for the RR states, one should generically expect 
two branches (rather than one) of gapless excitations in 
the ES, with different velocities. This is of course a well- 
known feature of the energy spectrum of a quantum Hall 
system with a physical edge. Actually, since we have 
shown that the real-space ES is the spectrum of a local 
field theory along the cut, most of the usual material on 
the edge spectrum can be safely used for the ES as well. 
This of course includes the case when there is an opera- 
tor that is less or equally irrelevant to the stress-tensor 
in the spectrum (for a recent discussion of the possible 
perturbations in the edge spectrum, see [50j). As an 
illustrative example that is not exactly a FQHE state, 
but that is closely related, one can look at the case of 
I ± i£-paired superfluids (spinless when I is odd). For 
those, the edge theory is known to correspond to ^ chiral 
Majorana modes at the boundary [56]. One generically 
expects to have a bilinear term that couples the different 
Majorana modes. Such a term in the bulk of the 1 + Id 
theory is relevant; it has scaling dimension h — \ < 2 
(this term is directly a bulk perturbation in the ES, we 
do not view it as a boundary perturbation along the edge 
of the sample). It leads to a splitting of the energy lev- 
els of Majorana modes [SOj. As should be expected from 
the generic arguments presented here, such a splitting 
should also appear in the real-space ES of these chiral 
superfluids. This was indeed found in Ref. [22) . 



3. RG analysis of the ES: the example of tfie Laughlin state 

Let us discuss the case of the Laughlin wavefunction 
in some more details now. First, we note that for dif- 



ferent particle number Ny^ ^ N'j^^ two states 



and 



always have a zero overlap. This rules 



out the possibility of having the operator e*"^'^)/^ or 
e-'y(^)/ ^ in the perturbation Sb{A) of the boundary 
action ( [4.15 ), or any other vertex operator, as they 
would allow non-zero overlaps between states with dif- 
ferent particle number. Thus, for the Laughlin wave- 
function, the most generic local perturbation along the 
boundary takes the form of a sum of polynomials in the 
(derivatives of the) U{1) current id(p{z). To avoid tech- 
nical issues caused by the extensive U{1) charge of the 
droplet, it is more convenient to use the shifted chi- 



+ i 



loe 



z. 



It is de- 
„ , where 



ral bosonic fleld ip{z) — ip{z) 
fined such that id(p{z) — \Jo + Sn^o 
J„ = 2^ / (""i-dfiC) is a Fourier mode of the origi- 
nal { i.e. not shifted) C/(l) current id(p{z) (see also sec- 
1), and only the zero mode is shifted Jq = 
/!/. This ensures that, acting on the CFT 



tion 
Jo - 



II B 



Nao 
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are normal-ordered. The coupling X^k} is of order 



l-lfelH hfcp 



vacumm with Nao + ANao charges, the eigenvalue of Jo 
is Jo \Nao + ANao) = ANao/V^\Nao + ANao), which 
is of order 1 in the scaling region, while the eigenvalue 
of Jo would rather be of order Nao- Leaving aside these 
technicalities for now, the most generic boundary pertur- 
bation has the form 

SbiA) = J2 hk} I \dz\ i^d'^^) . . . (z5'=-^)(z) 
{k} -^l-H^ 

(4.25) 

where the sum runs over the finite sets {k} = {ki , . . . , kp} 
with ki,...,kp > 1. The polynomials in id(p{z) 

{k} 

" ^ . The least irrelevant operator is actu- 
ally marginal: it is the current id(p{z) itself, which has a 
scaling dimension 1. The next contribution corresponds 
to the stress-tensor T{z) = {id(p)'^{z), with scaling di- 
mension 2. There is another operator with scaling di- 
mension 2: id'^(p{z). Then there are three operators with 
scaling dimension 3, five operators with dimension 4, and 
so on. The same discussion holds for Sh{B) in the part 
_B, and we finally see that the operator S'es itself must 
have the form ( 4.25[ ). 

Further restrictions can be put on the operators that 
can appear in Sb{A) and Sb{B), and S'es- Usually, in 
the literature, the cut between A and B is chosen such 
that the bipartition is symmetric with respect to the ex- 
change of A and B. This is the case in particular for 
the sphere, when A and B are chosen as the northern 
and southern hemispheres. For the Laughlin wavefunc- 
tion, the symmetry under the exchange of A and B im- 
plies a symmetry the transformation of the U{1) current 
id(p{z) I— >■ —id(p{z), which prevents all the terms that 
contain an odd number of factors (iOkjifiz)), such as 
idip{z), id^ipiz) or {id(pf'{z), from appearing in the op- 
erator S'es. Then the least irrelevant operator is indeed 
the stress-tensor, leading to a splitting of the pseudo- 
energies of order v/R ^ I/^/Naq, corresponding to the 
linear dispersion relation. The next leading corrections 
are due to higher order terms, all of which contain an 
even number of factors {id^^ip), such as {id^(f)^{z), or 
{id(p)'^{z), both with scaling dimension 4, leading to ad- 
ditional splitting between the pseudo- energies at order 
{^/Pq~^ / RY — There is no correction of or- 

der I/Nao, because the only possible operator of scaling 
dimension 3 that is even under id(p{z) i— > —id(p{z) is 
{idcp){id'^(p){z), which is a total derivative (of the stress- 
tensor). 



D. Particle partition entanglement spectrum in the 
scaling region 

As mentioned in the introduction of this paper (section 
ID I, the particle partition (PP) has also been considered 
in the literature. Formally, this bipartition is very close 
to the RSP, which allows to treat the ES for PP in a 
way that is almost identical to the one for RSP [TJl [H] . 



In particular, our definition of the "scaling region" (sec- 
tion IV B ) can be extended straightforwardly to PP. On 



the sphere, using the SO(3) rotational invariance, the 
Schmidt eigenvalues and eigenstates for PP can be clas- 
sified according to their total angular momentum AL^ 
(rather than the ALf component as we did for RSP). We 
also fix some finite ratio Naq/N (usually A^^o = N/2 for 
N even). The scaling region then corresponds to the set 
of all Schmidt eigenvalues/eigenstates which correspond 
to values of ANa and AL^ that remain bounded (< M 
for some fixed M) when one sends N (and Nao at the 
same time) to infinity. 

Although we just emphasized that the total angular 
momentum La is the most natural quantity to use to 
classify the Schmidt eigenvectors/eigenvalues, since those 
do not depend directly on Lf, it is still worth looking at 
what happens at large Lf. Indeed, the Schmidt eigen- 
state corresponding to the largest value of Lf is nothing 
but a state where Na = N/2 particles with "pseudospin" 
A fill the northern hemisphere, while the ones with "pseu- 
dospin" B fill the southern hemisphere. Therefore, PP at 
large values of Lf is very similar to RSP with a cut along 
the equator. We will use this below to analyze the scaling 
behaviour of the ES. After stereographic projection, the 
equator corresponds to a circle in the plane where the 
coordinates z live. We note R the radius of this circle. 



L Completeness relation, corrections to scaling and 
entanglement spectrum 



The decomposition ( |4.7| , obtained from a complete- 
ness relation in the CFT space, is still valid for PP. More 
precisely, for any fixed A^i, the groundstate wavefunction 
'i'{zi, . . . , zn) can be decomposed as 



N 



Na 

{NA,k\'[[a{z,)\0) 

3 = 1 



(4.26) 



where we have used the completeness relation Ijsf^ = 
J2k \NA,k) {NA,k\ in the Na subsector. This equation 
can be derived using the procedure of section |IV A| re- 
stricting the positions of the operators inside/outside a 
disc of some radius as in Fig. [sj Since the correlators 
{Na, k\ Uj aizj) |0) and (A^| JJr ai^,) \Na, k) in the right- 
hand side are all analytic, the identity remains true for 
any positions of the points ZiS and z^-'s, without the re- 
quirement that they lie inside/outside a disc of some ra- 
dius. Thus for PP, we have a decomposition identical to 

(4.7), where the states vPA^ 1.1// and 
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replaced by 



4- 



{NaM 



72.exp 



B 

\NA,k) 



PP 



. "AO 
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yPP 



{NA,k\ 
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A / e^^^'^^/^d^z a(z)® c^(z) 
1 



|0) |0)) 



X (4.27a) 
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{N\ 



7?.exp 



A / e^^^'^^/^d^^ a{z)(g)c\z) 



R 



(4.27b) 
|iV^,fc) |0)). 



Note that the integration domain in the exponential is 
C, rather than A or B for RSP. Again, the value of the 
coefficient A is non-universal, and must be tuned like for 
RSP. It depends on Nao (or R) and is not the same as 
for RSP, because of the integration domain which is C 
rather than A or B. The normalization factos for PP are 



^pp _ 



7PP 
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— / Y[e''^^-^U'z.\{NAo\llaiz.)\0)\ 
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(4.28a) 
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ZZo = ^ / n e^^^-^-^d'^. I {N\ n a{z,) 



AO) 



Again, note that the integration domain is ' 
A or B. 



(4.28b) 
rather than 



The states 



\XlA 

^{NaM 



and 



PP 



^\NA,k) 



SP 



are ex- 



actly the edge states constructed in section |li B| In the 
scaling region, they all correspond to edge excitations of 
the droplet of radius R. Then the entire discussion of sec- 
tions |IYbJ |IVC] and |IVC2] applies as well for PP in the 
scaling region. The conclusion is that the PP-ES is also 
the spectrum of a local operator along the "cut" , namely 
the circle of radius R. Up to a global constant shift, the 
pseudo-energies are the eigenvalues of , related to the 
boundary perturbation of the action of the C FT S^^ [A) 



and Sf^{B) by the same formula as (4.18). Then the 
determination of the PP-ES boils down to a discussion 
of the operators that can/cannot appear in the boundary 
perturbations S^^{A) and S^^{B), and in the operator 
, exactly like in the case of RSP. Once again, this 
needs to be done case-by-case. 



2. Absence of the stress-tensor for particle partition 

The most striking difference between PP and RSP is 
the following. For RSP we argued that the ES usually 
obeys a linear dispersion relation because of the presence 
of the stress-tensor of the chiral CFT T{z) (or possibly 
T'^W(^z) and T^{z) separately) in the set of perturbing 
operators along the cut. This is expected to be a very 
generic boundary perturbation, and as we highlighted, it 



is related to the extrapolation length v that appears very 
often in the field of surface critical phenomena. However, 
the stress-tensor cannot appear in the particle ES. One 
simple reason for this is that, on the sphere, the 5*0(3) 
rotational invariance implies that the pseudo-energies de- 
pend not directly on Ai^ (which would be related to the 
Virasoro generator Lq in the CFT) , but rather on the to- 
tal angular momentum /S.L^. The stress-tensor (as well 
as its derivatives) is not a perturbing operator for PP. 
Thus, when there are no operators that are less irrele- 
vant than the stress-tensor, one has to look for operators 
with scaling dimensions large than 2. The determination 
of such operators depends on the CFT that underlies the 
wavefunction. 



3. The example of the Laughlin state 

In the case of the Laughlin wavefunction, the analy- 
sis can be pushed further. As for RSP, the ope rator s 



S'PP(^), S^^{B), and S^^ have the generic form ( |4.25 l 
Again, some additional restrictions are imposed 
erators in S'£g have to be even under id(p{z) i— ^ - 



ae op- 
id(p{z). 

Since the stress-tensor and its derivatives are not here 
either, the least irrelevant allowed operator has scaling 
dimension 4, and it leads to a splitting of the pseudo- 
energies of order ^ /Rf ~ when N, Nao 
oo. Hence, this operator is (id"^ (p)"^ {z) . Before we ana- 
lyze the ES that follows from this operator, let us give 
a slightly more detailed argument which shows how this 
operator appears. 

For large Lf, the Na = N/2 + AN a particles fill 
the northern hemisphere, while the Nb = N/2 — ANa 
particles fill the southern hemisphere. The free ener- 
gies — log and — log Z^^^ are the ones of a one- 
component plasma (in its screening phase as long as l/v 
is not too large) on the sphere, with a boundary along 
the equator. Generically, these two free energies should 
both contain a surface tension contribution, proportional 
to the perimeter of the cap that they fill. More pre- 
cisely, such a term does not appear at v = 1 (this is 
of course consistent with the fact that all the pseudo- 
energies are degenerate in that case), but when v > \ 
it is known that the surface tension is not zero [571 ISH] ■ 
Now we imagine that the droplet is slightly deformed, 
namely that its boundary is a curve parametrized by 
= I -\- 59{(j)), for a small displacement 59. Here {9,<j)) 
are the polar coordinates on the unit sphere. Then the 

length of the interface is d(t)\J &\v? 6 + {d9/d(j)Y — 

27r + i / {{d69/d(j)f - (60)^). The displacement of 
the interface is proportional to the current id^Lp(4>) (see 
[18j). Thus, we find that the surface tension term is 
given by the integral along the equator of the operator 
{id'^(fi)'^ {(j)) — {id^0)^ {(j)) . One can check that this is the 
same operator as {id'^(pf'{z) when one maps the sphere 
back to the plane, taking into account the Jacobian of 
this transformation carefully. Looking at the surface ten- 
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sion term on the sphere rather than on the plane has 
the advantage of making the 5*0(3) rotational invariance 
more transparent. This term appears both in Sf^{A) 
and S^^{B), with no obvious cancellation between the 
contributions coming from A and B, therefore it must 
also appear as well in S^g . We analyze the spectrum of 
the operator that we just found next. 

When it is written in terms of the modes of the 
(shifted) U{1) current idif — KJq + X^n^o 
zero Fourier mode of {id'^(pY{z) leads to 



cPP _ 
'-'es — 



C 



N' 



3/2 
AO 



fe>l 



(fc^ — l)J-kJk 



Oil/N 



lo) 



(4.29) 

where C is some undetermined positive constant of order 
0(1). This constant vanishes for the integer quantum 
Hall effect (Laughlin with i/ = 1), but it is non-zero 
as soon as < 1. As explained previously, the par- 
ticle ES is nothing but the spectrum of Sgg, up to a 
global additive shift that ensures the normalization of 
the reduced density matrix trp^ = 1. The spectrum 
of the operator (4.29) has several interesting features. 



The eigenstates of S'es at the order l/N'^^ are all the 
states nfc>i •^-k \^Ao + ATV^a), for any finite set of posi- 
tive occupation numbers n^'s. Then, for each finite set of 
non-zero occupation numbers, the splitting between the 
corresponding pseudo-energy and the ground state is, at 
the leading order 



the trend ends at this k; the next one fc = 7 is degener- 
ate with that at fc = 6. These levels (with fc < 6) should 
agree with setting nk — I, and others zero, for fc = 1, 
2, . . . , 6, in the above formula, and so lie on the cubic 
curve = {C/N%^){k^ - fc) plotted in the Figure. We 
find a reasonable agreement, as shown, for C — 0.5144 
at this size. At the same time, the lowest pseudoen- 
ergies for each Na should lie on an inverted parabola, 
and the coefficient should be the same according to the 
above formula (4.301. A plot is shown in the inset on the 



Figure, with the same parameter value, and the agree- 
ment is reasonable. If the two curves are fitted to these 
forms independently, with one parameter for each, then 
the best parameter values are found to be within about 
10% of each other. Further features of the spectra also 
agree with the above form. Let us note first that the 
cutoff at fc = Na noticed above makes sense as the stan- 
dard cutoff from finite size effects in the edge spectra; the 
sums of powers X^ilfi j which create these elementary 
excitations if trial wavefunctions are used, become alge- 
braically dependent on the lower ones if fc > Na- Levels 
in the spectra other than the extreme ones fitted above 
can be roughly explained as sums of multiple excitations, 
in rough agreement with the formula. In particular, the 
banded structure that is apparent in the Figure emanat- 
ing from each of those extreme levels can be understood 
in each case as excitations with lower fc' < fc added to 
the fcth extreme one. Similar structures and trends are 
seen in smaller sizes, and for Na < N/2 (not shown). 



AC = - 
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N 



3/2 
AO 



{ANa)' 
2u 



-Y,nk{e-i)k 



fc>0 



(4.30) 



This spectrum possesses a few noticeable features. First, 
looking at the lowest pseudo-energy in each AN a sec- 
tor, which always corresponds to no non-zero occupa- 
tion number rifc, one gets the inverted parabola observed 



numerically in |15j : A^ = — 



C {ANaV 

3/2 



N 



2v 



Second, for 



fixed AA^^, the spectrum of ^es is in general highly 
degenerate. This is due to the fact that A^ does not 



depend on ni 
all the states 



Thus, in the Schmidt decomposition. 



\llA 



\XlA 



PP 

contribute 



\ijA 



with 



PP 
the 



PP 



\ijA 

(A^|.'f//pp- 

same pseudo- 



energy (at least at this order), leading to a large de- 
generacy in the ES. Such a large degeneracy is to be 
expected, since the 5*0(3) rotational implies that the 
pseudo-energies appear in multiplets. 

In Fig. § we compare the formula (|4.30|) with nu- 



merical results for the particle ES of the Laughlin state 
with = 1/3. We plot low- lying pseudoenergies ^ versus 
for the scaling region near L^g, for Na = N/2 and 
A^ = 12 particles. The range of pseudoenergies included 
is a little larger than in Ref. [15], to show the structure. 
It is interesting to look at the upper envelope of the levels, 
namely the maximum ^ for = Lf^ — fc for each fc = 1 , 
2, .... For fc < 6 = Na, those levels fall on a curve, but 
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FIG. 6. Plots of PP pseudoenergies ^ versus Lf for Na = N/2 
and versus Na for the N = 12, v = 1/3 Laughlin state. Here 
A*,^ — ^ — Nln2. Main figure: levels in the scaling region ver- 
sus Lf. The dotted curve is a one-parameter fit of a cubic to 
seven levels, explained in the text. Inset: lowest pseudoenergy 
for each Na; the values for Na ~ 0, 12, which are ^ = A''ln2, 
are omitted. The curve is an inverted parabola, with the same 
parameter value as the main figure, as explained in the text. 



Let us finally mention that, interestingly, our results 
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for the Laughlin wavefunction can also be related to a 
series of papers by Zabrodin and Wiegmann, who have 
studied the behaviour of the droplet as the potential 
V{z, z) is varied, using the so-called "loop equation" 
[IHllMlllSl- The details of the calculation that lead to the 
operator (id'^(pf' using this alternative approach are be- 
yond the scope of our paper. Here we just want to point 
out that one should be careful when trying to relate their 
results to ours. Indeed, the fact that the surface tension 
for the one-component plasma vanishes ai v = \ has been 
the source of some confusion. In particular, the term cor- 
responding to the surface tension, which is non-zero when 
I' < 1, is actually missing in the result stated by Zabrodin 
and Wiegmann in |48j, although they do consider it in 
their calculations (it appears as a zero mode of the loop 
equation in their paper). This term, however, should be 
there, as pointed out for instance in [57l[58]. We note 
that, in the FQHE literature, the possibility of a surface 
tension term was also quickly discussed in the appendix 
ofRef. EI]. 



V. DISCUSSION: ENTANGLEMENT, TRIAL 
WAVEFUNCTIONS, CONFORMAL BLOCKS 
AND TENSOR PRODUCT STATES 

Since the seminal work of Laughlin [2] , it is known that 
an efficient approach to the theory of the FQHE con- 
sists in making various proposals for trial wavefunctions, 
which are not directly related to the eigenstates of the 
original Hamiltonian, but which still accurately capture 
some of its physical features. This approach is of course 
not specific to the FQHE effect, and other strongly cor- 
related quantum systems can be studied with variational 
wavefunctions. Some powerful techniques have been de- 
veloped for onc-dimcnsional systems, especially for spin 
chains, which focus on some class of trial wavefunctions 
called Matrix Product States (MPS) [.62 . As a simple 
example of these, let us consider an open spin-1/2 chain 
with N sites. An MPS for this system can be constructed 
as 

m = (5.1) 

{Vont\ A{sn) ■ ■ ■ A{s2)A{si) \vin) \ S1S2 ■ ■ ■ Sn)) 

Si,....SN=t,i 

where the states | S1S2 . . . sn)) with Si =t,i span the 2^- 
dimensional Hilbert space H = {t,!}*^^, and A{t)Mi) 
are both D x D matrices acting in an auxiliary space Vd- 
The two states \vin) and \vout) belong to the space Vd- 
The choice of these two states and of the matrix A fixes 
the state | 5')). 

An MPS possesses peculiar entanglement properties: 
dividing the system into two subchains of length Na and 
Nb (with Hilbert spaces -Ra = {t,i}®^-* and T^b = {t 
,1}^^'^), one can easily check that the Schmidt decom- 
position corresponding to the bipartition H = Ha 'S)'Hb 



is of rank D at most. This simply follows from insert- 
ing a completeness relation in the space Vd (namely 

Ivd = J2k=i \^k) {vk\ with {vk\vk') = SkM') between the 
matrices A{sn^^i) an A{sna)- This is a striking feature 
of an MPS, as it is radically different from generic states 
in H which would typically have Schmidt decompositions 
of rank min(dim'H^, dimHs). In that sense, the set of 
MPS in H for some fixed D can be thought of as being 
a submanifold (of viewed as a manifold) of minimally 
entangled states. In particular, when D — 1, the space of 
MPS is nothing but the submanifold of separable states. 

Starting from a Hamiltonian H for our spin-1/2 
chain, one should minimize the expectation value 
((*! H l^f)) /((*|*)) within the set of D-dimensional 
MPS, which would give some approximation to the true 
ground state of H. By increasing the dimension D of 
the auxiliary space Vd, one can approximate the true 
ground state of the system with arbitrary accuracy. MPS 
are at the heart of numerical methods such as the Den- 
sity Matrix Renormalization Group (DMRG) [63], and 
have also been used recently as a key theoretical tool in 
the successful efforts made by several groups to classify 
one-dimensional interacting topological phases of matter 
[Bm67j . Their higher-dimensional analogs, the TPS (also 
called "PEPS" for Projected Entangled-Pair States [55] ) 
are currently being used in attempts to classify higher- 
dimensional topological phases (see for instance [69]). 
The basic idea of such a program is that any gapped 
phase of matter is believed to be accurately represented 
by one (or more) trial state, like an MPS for D large 
enough, or some higher-dimensional generalization (see 
for instance [70]). Such a trial state is not the ground 
state of the actual physical Hamiltonian, but it is the 
ground state of some local Hamiltonian [3 [7T1 [Z2] which 
is believed to be adiabatically connected to the physi- 
cal one. Then the topological properties of the phase, 
namely the properties which are robust as one varies the 
different physical parameters within the phase are cap- 
tured by the trial state, and the trial state has the enor- 
mous advantage of having some nice mathematical struc- 
ture, which makes analytic and/or numeric calculations 
tractable, even for large system sizes. 

We would like to point out that the states given by the 
MR construction have a structure that is very similar 
to the one of an MPS or a TPS. Essentially, it can be 
viewed either as an MPS with D = 00, or equivalently 
as continuous version of a TPS. Indeed, we note that the 



conformal block in (1.6 1 is nothing but 



(fout|a(2A')...a(zi)|win) 



(5.3) 



where |win) — |0) is the vacuum of the CFT, and \vout} = 
\N) is the vacuum with N charges. The structure of 



this correlator is clearly analogous to the MPS (5.11, 



m = Y,e-^^/'\^A,))<E>\^B,)) 



(5.2) 



with a space Vd which is now the (infinite-dimensional) 
Hilbert space of the CFT. The use of conformal blocks as 
infinite-dimensional MPS has recently been advocated in 
Ref. [73] ■ In particular, in this work, conformal blocks 
are used as trial wavefunctions for critical (gapless) spin 
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chains. We note that, in the context of the FQHE, such 
wavefunctions are nothing but the ones given by the MR 
construction (section I C I . The fact that Vjj is infinite- 
dimensional does not mean than the particular entan- 
glement properties emphasized above completely breaks 
down, since in each given angular momentum subsector 
(that is, in each Lq eigenspace) the auxiliary space is fi- 
nite, and thus the Schmidt rank in that sector is finite. 
This is exactly the property that was observed by LH, 
namely that particular trial wavefunctions for the FQHE 
have an infinite entanglement gap. We emphasized the 
connection between the so-called "level-counting" in each 
angular momentum sector, the conformal block struc- 
ture, and the completeness relation in the auxiliary space 
in section HV Al 

Another way of looking at the trial wavefunctions for 
the FQHE is to see them as continuous TPS in 2d. Let 
us imagine that the ground state in a two-dimensional 
gapped phase of matter can be approximated by some 
TPS. For a continuous system of identical itinerant par- 
ticles, a natural way of searching for a TPS approxima- 
tion of the ground state would be to discretize the space, 
namely to put the particles on a lattice. For simplicity, 
let us assume that the particles are identical hard-core 
bosons. Then on each site i of the lattice, either there is 
a particle, or there is none. There is one tensor associated 
to each these two situations, Al „ and Bl „ , with 
p indices corresponding to the p edges that end at the 
site i. Each of the ei can take values from 1 to D, where 
D is the dimension of the auxiliary space Vd- Then, 
in the ground state, the amplitude of the configuration 
I 0)) with n bosons on the sites zi, . . . , zat (for N 



■ ■ c. 



sites in total) is given by the contraction of all the edge 
indices of a product of tensors A {n times) and B {N — n 
times). Then we sum over all the possible numbers of 
particles n, and we get the trial wavefunction: 



I*)) = try. 



n ( 



A'(g}l + B' 



|0)) (5.4) 



where the trace is taken over the auxiliary space Vjj 
for each edge of the lattice, and we have omitted all 
the edge indices. For each fixed configuration of edge 
indices {ei,e2,...}, one may interpret the amplitude 
A^A"^ . . . A^ (where we have dropped the edge indices) 
as an amplitude. The same thing can be done with in- 
sertions of tensors i?*'s in some places on the lattice. It 
is now clear what the natural extension of a TPS should 
be when we replace the lattice by the continuum (here 
the two-dimensional plane, or sphere). The set of edge 
indices {ei,e2, . . . } should be replaced by local degrees 
of freedom, which we note h(x). A configuration of de- 
grees of freedom has an amplitude Yii • Iii the con- 
tinuum, the product on each point "Ox" more con- 
veniently written as the exponential of a sum, which 
becomes an integral. Thus, to a given configuration of 
variables {/i(x)}, we associate the amplitude e-l' ''■^^^'^'^^^^ . 
The continuous analog of the tensor i?*, which modifies 



the amplitude AK . . A\ . . A'^ to AK . . B\ . . A^ , is an 
operator 0[h{x.)]. Then the continuous TPS takes the 
familiar form 

I*)) = J [d/j]e^'^'^^['*('')]e°[''('')l®='('')|0)), (5.5) 

up to normalization. In other words, the am- 
plitude of the configuration c'^ (xi) . . . c'''(x„) 1 0)) is 
/[rf;i]e/'''''-^[^WlC>[/i(xi)] . ..0[h{^n)]- This, of course, 
is nothing but the correlator of a set of local operators in 
a local field theory. 

We thus conclude that, searching for a TPS-like Ansatz 
for a continuous system of mobile identical particles, one 
is naturally lead to considering local field theories, and 
correlators of local observables in those. This remark, of 
course, can be extanded in a straightforward way to the 
case of particles with additional degrees of freedom (such 
as spin), by using one operator O per degree of freedom 
(for instance O-f and Oi, with an appropriate SU{2) sym- 
metry, for spin 1/2). Although lattice TPS may appear 
as more useful in numerical simulations, working with 
their continuous analogs can be helpful in a few cases. In 
the FQHE in particular, and more generally in phases of 
matter with breaking of parity and time-reversal symme- 
try, one may encounter various issues by trying to search 
for lattice TPS. One reason for that is the following. For 
fermionic systems, the field theory/lattice model that one 
uses to contruct the trial wavefunction must be fermionic 
as well. If one wants a trial state which breaks parity and 
time-reversal symmetry, one is naturally lead to consider 
chiral fermions. It is a notoriously difficult problem to 
construct lattice models which give rise to chiral fermions 
|74j . Working directly in the continuum allows to circum- 
vent such problems. There are plenty of well-known ex- 
amples of field theories of chiral fermions. Such examples 
include, of course, some chiral CFTs. The use of confor- 
mal blocks as trial wavefunctions is thus natural in that 
context. The trial states for the FQHE that are given by 
conformal blocks may all be viewed as continuous TPS. 

To conclude this section, we note that the results given 
in this paper may be reformulated more generally for 
the larger class of (lattice or continuous) TPS. For all 
of these states, one can construct some wavefunctions for 
the edge excitations in a way similar to ours. These states 
are also the ones that appear in the Schmidt decompo- 
sition of the ground state. In general, two important 
questions are (i) how to compute the overlaps between 
these "edge states", and (ii) whether the entanglement 
Hamiltonian is isospectral to a local pseudo-Hamiltonian 
along the cut (up to distinctions of a finite number of 
topological sectors, which by definition are not local). A 
natural conjecture is that the property (ii) holds as soon 
as the connected correlations of local operators are all 
short-range in the bulk. In the TPS context, we note 
the existence of a recent work along these lines (Ref. 
[75]), which provides numerical evidence for the local- 
ity of the pseudo-Hamiltonian in a number of specific 
two-dimensional lattice models. 
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VI. CONCLUSION 

In summary, we have studied the trial states for the 
FQHE that are given by conformal blocks. We have con- 
structed the wavefunctions for to the gapless edge ex- 
citations that are associated to those. We define those 
directly, relying on the conformal block structure: to each 
CFT state we have attached a physical states | 
We did not address (physical) special Hamiltonians, in 
contrast with most previous approaches. In all known 
cases, however, and in particular for Laughlin and the 
MR Pfaffian state, these edge states agree with the ones 
constructed previously, for instance in [ISJ . 

Then, we have studied the overlaps between these edge 
states, relying heavily on the assumption that all the con- 
nected correlations of local observables are short-range in 
the bulk, a property dubbed generalized screening after 
[53] . The techniques involved in the study of the overlaps 
were the ones of boundary CFT [331 SO] • We used a gen- 
eralized plasma analogy to compute these overlaps. The 
basic idea can be summarized as follows: the chiral CFT 
that underlies the trial wavefunctions lives in the complex 
plane C. The complex conjugate of the wavefunctions are 
obtained from an anti-chiral copy CFT. In the thermo- 
dynamic limit, there is a simply connected region Q G C 
filled by the N particles: the droplet. Inside the droplet, 
generalized screening means that the non-chiral theory 
CFT (g) CFT is perturbed by a local operator a{z)(^a{z), 
and flows towards a massive fixed point under RG. Then 
we are left with a massive field theory inside the region Q, 
and a massless theory CFT® CFT in C\ri. At the bound- 
ary of Q, the massless theory obeys a conformal boundary 
condition. We have determined this conformal boundary 
condition (a(z) oc a^{z) at the boundary), and shown 
how it implies that the (antilinear) mapping from the 
CFT space to the physical space \v) | becomes 
an isometric isomorphism in the thermodynamic limit 
N ^ oo. This is a precise bulk/boundary correspon- 
dence, for which we provided some elementary numerical 
checks (Monte-Carlo) in the case of the Laughlin and the 
Moore-Read (Pfaffian) states. In particular, one conse- 
quence of this isometric isomorphism between the CFT 
space and the physical space of edge excitations is that 
the CFT must be unitary. We have also shown that an- 
other aspect of this bulk/boundary correspondence holds 
as a simple consequence of this boundary CFT formal- 
ism: the equal-time correlators (evaluated in the ground 
state) of local operators along the edge are equal to the 
analogous correlators in the chiral CFT. 

We then tackled the entanglement spectrum of the 
ground state for real-space partition. A key point was 
the use of a completeness relation in the CFT space, that 
leads to a natural decomposition of the ground state | 5')) 
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In the 



as a sum of terms of the form 

limit Na,Nb — >, because of the isometric isomorphism, 
all these states are orthogonal, and therefore this natural 
decomposition is really a Schmidt decomposition. All the 



pseudo-energies are degenerate in that limit. More seri- 
ously, most of the interesting features of the ES are hid- 
den in the subleading corrections. We sketched how such 
corrections to scaling can be naturally apprehended with 
renormalization group group arguments: in reality, the 
massless theory CFT (g) CFT that lives in C\r2 is subject 
to a perturbed conformal boundary condition along the 
edge. Therefore, one needs to discuss what perturbations 
are allowed by the different symmetries to be present at 
the boundary. A careful analysis of these perturbations 
leads to a deformation of our isometric isomorphism, due 
to the corrections to scaling. Of course, the latter need 
to be taken into account in the orthogonalization of the 
states that appear in the Schmidt decomposition. These 
corrections to scaling are thus reflected in the entangle- 
ment spectrum. Our most significant result here is that 
the ES, as a consequence of generalized screening, is the 
spectrum of a sum of (integrals of) local operators along 
the cut between A and B. We also showed that the ES 
for particle partition (PP) can be analyzed in a similar 
fashion. 

Finally, we have emphasized that our results suggest 
relations between the formalism of wavefunctions given 
by conformal blocks and Matrix Product States (in Id) or 
Tensor Product States (in higher dimension) , noted MPS 
or TPS. We have argued that the continuous analog of a 
TPS (defined on a lattice) is nothing but the correlator of 
local operators in a local field theory. Working in the con- 
tinuum rather on the lattice has a number of advantages, 
in particular when one is dealing with the ground state 
of a phase of matter with broken inversion/time-reversal 
symmetry, such as the FQHE. In that case, assuming 
that good trial wavefunctions can be built as correla- 
tors in a local field theory, one is naturally lead to look 
for chiral field theories, and correlators of local fields in 
those. A good illustration of this is the case of com- 
plex paired superfluids, for which the BCS ground state 
can always be viewed as a correlator in some chiral free 
fermion theory. Looking for a lattice TPS approximation 
of the ground state is then equivalent to constructing a 
lattice discretization of a chiral field theory. Then one 
runs into trouble, as various problems arise, in partic- 
ular the fermion doubling problem |74j . Starting from 
the continuum, namely from a field theory, to construct 
a trial wavefunction, is a natural way of circumventing 
such issues. This justifies, at least heuristically, the ap- 
proach pioneered by MR, namely the use of conformal 
blocks as trial wavefunctions for the FQHE. 
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Appendix A: The topological term 7 appears in all 
the pseudo-energies 

In this appendix we give one argument to support the 
fact that the topological entanglement entropy 7 appears 
in all the pseudo-energies in the scaling region (see sec- 
tion IV B I . It is adapted from the argument sketched by 



Kitaev and Preskill at the end of Ref. 

We have shown in this paper that the ES (for RSP) 
is, up to a constant shift, given by the spectrum of 
^Lq + o{y/])o~^ / R). The constant shift is fixed by the 
requirement that trp^ = 1. To compute this shift, let 
us start by dropping the terms of order o{y/pQ~^ / R). 
Within this approximation, the reduced density matrix 
PA is isospectral to 



1 



(Al) 



acting on a module V over the chiral algebra A; usually 
the module of the identity. The normalization factor is 



Z = gf24 X try q 



(A2) 



with q = exp {—'^'^2nR}- "^^^ central charge of the CFT 
is noted c. The object in the right-hand side is called a 
character in the CFT literature p7H29] . 

The constant shift we are looking for is the order 0(1) 
piece in the expansion of logZ = a2TTR — j + o{l / R) . We 
thus need to obtain the asymptotic behaviour of log 2 
when i? ^ u ~ v^Po~^, namely when q — > 1. At first 
sight, this may not look obvious, because the sum di- 
verges in that limit, however it is a very standard ex- 
ercise in CFT. The trick is to use the modular proper- 
ties of the characters, namely to express the character 
x{q) = trv'(7^°~5i as a function of g = e^"^^^^, and 
then take the limit g — >■ 0, which is well-defined. In gen- 
eral, the characters obey the following rule 



Xaiq) = J2^^ ^biq) 



(A3) 



where a and b label the different irreducible modules over 
A, and is the so-called modular S-matrix (note that 
this is not the same object as the shift operator S is 
In the case when V is the module of the 



section 



IIIB) 



identity, the leading contribution comes from the char- 
acter of the identity module itself, noted xi (and the 



corresponding diagonal element of the modular 5-matrix 
noted Si) 



logXi(9) log 



log xiiq)] 



24 



log q + log Sl- 



(A4) 

(A5) 
(A6) 



Since logZ = logxi((7) + ^ logg, we find in the end 



logZ = a2TrR- -f + o{v/R) 



(A7) 



with 7 = — logiS^ as expected. The coefficient a is not 
universal, as it involves the inverse length 1/v. The ar- 
gument given here, based on modular manipulations, can 
be easily extended to the case when anyonic excitations 
are present in parts A and B. More details can be found 
in ITO1I76]. 

We do not address the robustness of 7 here. Since 
it is a dimensioneless constant, it cannot be affected by 
irrelevant perturbations of the CFT. One can make some 
sanity checks in a few cases, for instance for a CFT with 
both a C/(l) sector and a statistics sector, which come 
with different velocities w;7(2) and v^. One finds that 7 
is invariant indeed. 

In conclusion, the topological term 7 must appear as 
the order 0(1) piece in all the low-lying pseudo- energies 
in the ES. Because of this, it also appears in the von 
Neumann entanglement entropy, but this is not specific 
to this entropy. It is a property of all the pseudo-energies 
independently, rather than a property specific to the von 
Neumann entropy only. In particular, one consequence 
of the presence of the term 7 in all the low-lying pseudo- 
energies is that all the Renyi entropies contain this term 
7 as well. This fact was noted already in [77] . 



Appendix B: The exactly solvable case of px ± ipy 
paired superfluids 

In this appendix we discuss the case of a p^; — ipy su- 
perfiuid treated in BCS theory [56], where the overlaps 
between the edge states and the real-space entanglement 
spectrum can be computed exactly. We reformulate some 
of the results of \22i in the language of the present paper, 
to illustrate the ideas developed in sections IC III and 

El 

For convenience, we consider the system on the cylin- 
der C. The BCS ground state of spinless particles is ex- 
pressed in terms of the pairing function g{wi,W2) 



1 



z 



(Bl) 



d^wi I d^W2 g{wi,W2) c\wi)c\w2) ] |0)) 



exp 



where w ^ x + iy is the complex coordinate on the cylin- 
der (y and y + L are identified). For a weakly paired 
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Px ~ iPy superfluid, the pairing function expressed in mo- 
mentum space behaves as 



g(k) oc 



1 



ik. 



(B2) 



V 



as |k| — ^ (see [53]). Coming back to real space, on the 
cyhnder C, the pairing function is 



Ij/2tt 



sinh 



27r(wi—W2) 



(B3) 



where /i is a complex parameter with the dimension of an 
inverse length, which is related to the density of particles. 
Without loss of generality, one can assume that fj, is real 
and positive. For this particular choice of the pairing 
function, the BCS ground state (Bl I can be rewritten as 



Note that the chiral correlator is still evaluated on the 
full cylinder, but the particles fill the left half-cylinder 
only. The normalization factor then corresponds to the 
partition function of the Majorana fermion field on the 
cylinder, with a mass term which is switched on in the 
left side only 



Si 



1 

2n 



2n 



d w tptp 



(B8) 

From now on we make use of the translation invariance of 
the system in the y-direction, by using a fcy-momentum 
basis. On the cylinder the field 'ip{w) can be expanded in 
Fourier modes 



tp{x,y) 



^^e'=«(-+'^)V-.„ (B9) 



1 



^ji:^ I d'w 'ilJ[w)^c'[w) I ) |0)) 

'(B4) 

where 'ip{w) is a free Majorana field on the cylinder with 
propagator (-0(wi)^/;(-u;2))c = l/(^ sinh ^^t"'^-'"'') ). 
The chiral correlators {^^(wi) . . .'iP{wn))c evaluated 
using Wick's theorem. Introducing an anti-chiral copy 
ip{w) of the Majorana fermion field, the normalization 
factor Z can be expressed as a correlator in the non-chiral 
theory Majorana Majorana 



7 = \ " V27r/ 

^ (n!)2 

n>0 ^ ' 



l\^{wi) . . . tp{Wn)lp{w[) . . . tp(w'„)) 

X ((c(z„)...c(zi)ct(zi)...ct(z;))) 

- J d~w 'ip{w)ip{w) 



(B5) 



The factor i in the last exponential appears because of 
the convention that the holomorphic/anti-holomorphic 
modes to anticommute {V'lV'} = 0- We also use the 
fact that only the correlators with an even number of 
fields ipi^i) ^re non-zero. The term in the last exponen- 
tial should be viewed as a mass term in the Majorana 
fermion theory, with the standard quadratic euclidean 
action 



S = 



1 

2^ 



tfw [ipdwi^ + i^dyj^p + ifj, ^ip] . (B6) 



a. Edge states 

We now focus on the Px + ipy superfluid filling only 
the left half-cylinder Ci corresponding to a: < 0. A trial 
wavefunction for this system is obtained by restricting 
the integration in the expression (Bl ) to the domain C/, 



1 



|0)). (B7) 



where the sum runs over all the admissible momenta, 
namely ky e ^ (Z + ^ ) for anti-periodic boundary con- 
ditions. The Fourier modes obey the canonical commu- 
tation relations {'0fc„i'0fcj,} — ^ky+k'^fl- The ground state 
trial wavefunction (B7| is then 



I T exp 



1 



(BIO) 



dx e'^y^ ij.ky^cUx) |0) |0)) 



where the symbol T denotes time-ordering for the imag- 
inary time X. The state |0) is the vacuum of the CFT, 
annihilated by all the mode V'fc„ with ky > (note the 
difference with the physical fermionic vacuum |0)) an- 
nihilated by the c(w)'s). The time-ordering operator T 
implies that the Cooper pairs of particles with momenta 
+ky and —ky must appear in a specific order, as discussed 
in 122 . 



As explained in section II B one can write down trial 



states for the edge excitations by replacing the out 
vacuum (0| by an excited state in the CFT, such as 
(0| 'ipkt'^k2 for some ki,k2 > 



1 



(Bll) 



(0| VfciV-fc. r e^^-y ^""^ ^--y^-lyi-) |o) 1 0)) 

The study of the overlaps between these states is equiv- 
alent to the one of correlators of of the fields ij) £^nd -ip 
in the right half-cylinder Cr {x > 0), in the presence of 
the massive phase in the left half-cylinder C; (x < 0). In 
the limit /i — > oo, one is left with a conformal boundary 
condition at .t = for the massless Majorana fields in Cr 



^(0,2/) = V(0,2/) 



(B12) 



which gives, in terms of the Fourier modes acting on the 
boundary state \B) 



^fe. \B) = ^-k,. \B) 



(B13) 
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The boundary state \B) can easily be written in term of 
the fermion modes 



\B) 



|0)|0). 



(B14) 



As explained in section pIID[ this implies that the over- 

(B15) 



laps in the fi ~> oo limit are given by 



(0|^,,...^, 

(o| tZ-fci . . . (o| ■ ■ ■ ipk- 
3| V'fel • ■■i'k'^i^k^ ■■■■>Pki |o) . 



\B) 



We have arrived at this result thanks to the fact that the 
massive phase of the theory Majorana® Major ana in the 
domain Ci can be reformulated as the boundary condition 
(B12). This is a particular case of the general argument 



of section HID which leads to the result (3.231. 



Since the + ipy-paired superfluid wavefunction (Bl 



is Gaussian, one can actually compute the overlaps (B15 ) 



exactly for any fj,, not only in the limit — >■ oo. This can 
be done as follows. For some given ky,k'y > 0, we have 



fi / dx e^y'^ Cky{x) 

J —oo 





dx' e^'y-' cl,ix')\^,Ci)) 



^ "ky.k'y 



where 



St,{ky) = log + {ky/nf + ky/f?j 



(B16) 



(B17) 



Here we have used a result of Ref. [22] (see formula (13) 
and (14) in that reference). The case with more fermionic 
excitations follows simply from Wick's theorem: in gen- 
eral, the overlaps at finite fj, are 



{0\4,k^...^kp 



= (0| Vfe; . . . ^k'^e-^'^^'^k, . . . V'fci |0) . (B18) 
where we have defined the boundary perturbation as 

Sbifi) - SM^kyi^-ky (B19) 

ky>0 

Of course, in the limit fi — >■ oo, the boundary perturba- 
tion goes to zero, and we recover the formula (B16|. 



b. Real-space entanglement spectrum 

The real-space entanglement spectrum on the cylinder 
(where A is the left half-cylinder Ci and B is the right 
half-cylinder Cr) can also be computed exactly, as pointed 



out in [55]. Because the left and right half-cylinders are 
symmetric here (such a symmetry was not assumed in 
Part IV I , the boundary perturbations for part A and part 



B are both equal to Sb{A) — Sb{B) — 5*6 (/i) defined in 
(B19). Following the discussion in section IV and in 
particular the formula ( 4.18[ ), we define 



(B20) 



and the entanglement spectrum is directly given (up to 
an additive constant shift) by the spectrum of the op- 
erator Ses = '^Sh{p), which is a free fermion spectrum, 
generated by the set of single-particle pseudo-energy lev- 
els 2S^{ky) for the admissible ky > 0. The constant shift 
can be computed easily using this fact. For example, if 
we focus on the smallest pseudo-energy ^q, corresponding 
to the cut ground-states in A and B (that is, without any 
insertion of a ipky mode) , then we have the exact formula 



-&/2 



1 



ky>0 



-2Sj_t {ky) 



(B21) 



where the product runs over all the admissible ky > 0. 
Then is given by the Euler-MacLaurin formula as L — > 
oo 

Co = a i - 7 + • ■ • (B22) 

where the coefficient a is the integral dky2\og{l + 
g-25,.(fc„)-j^ Here the topological term is actually 7 = 0, 
in apparent contradiction with the claim made in sec- 



tion IV that this term is the topological entanglement 
entropy. However, the topological entanglement entropy 
itself can be computed exactly by the same method, us- 
ing a similar Euler-MacLaurin expansion, and the result 
would also be 0, again in contradiction with the expected 
result 7 = log 2. This comes from the fact that we have 
not restricted the parity of the particle number in part 
A (or equivalently B). The topological term can only 
appear in that case, and it does if we restrict the number 
of particles in A to be, say, even (see also the discussion 
in [H]). Thus the result is indeed in agreement with the 



general claim in section IV The set of pseudo-energies is 
then given by 



^ki,k2 



P 

E 



2S^{kj) 



(B23) 



where p is assumed to be even, and ^ L — log 2 + . . . . 
The same result holds if we restrict the particle number 
in A to be odd. 

Finally, to conclude this appendix, we note that, by ex- 



panding the formula (B19) in powers of fcj,//i, and going 



back to real space, the operator S^^ is indeed the integral 
along the cut between A and B of some local operators 
in the CFT. As expected from the discussion in section 
JVj the first term in ^es is nothing but the stress-tensor 
itself 



•Ses — 



27r 



my) 



(B24) 
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and the length l//x, which is the only microscopic length extrapolation length discussed in the main text (up to 
scale in this model, appears as a coefhcient. It is the some numerical factor). 
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